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Abstract 

We study the scaling limits of three different aggregation models 
on l/: internal DLA, in which particles perform random walks until 
reaching an unoccupied site; the rotor-router model, in which particles 
perform deterministic analogues of random walks; and the divisible 
sandpilc, in which each site distributes its excess mass equally among 
its neighbors. As the lattice spacing tends to zero, all three models 
are found to have the same scaling limit, which we describe as the 
solution to a certain PDE free boundary problem in . In particular, 
internal DLA has a deterministic scaling limit. We find that the scaling 
limits are quadrature domains, which have arisen independently in 
many fields such as potential theory and fiuid dynamics. Our results 
apply both to the case of multiple point sources and to the Diaconis- 
Fulton smash sum of domains. 
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1 Introduction 

Given finite sets A,Bc Z'^, Diaconis and Fulton [7] defined the smash 
sum A (B B as a certain random set whose cardinality is the sum of the 
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Figure 1: Smash sum of two squares overlapping in a smaller square, for 
internal DLA (top left), the rotor-router model (top right), and the divisible 
sandpile. 



cardinalities of A and B. Write A Ci B = {xi, . . . ,Xk}. To construct the 
smash sum, begin with the union Cq = AU B and for each j = 1, . . . , fc let 

Cj = u 

where yj is the endpoint of a simple random walk started at Xj and stopped 
on exiting Cj-i. Then define A(BB = Ck- The key observation of [7] is that 
the law of A(BB does not depend on the ordering of the points xj . The sum 
of two squares in overlapping in a smaller square is pictured in Figure [ij 



In Theorem 1.3 below, we prove that as the lattice spacing goes to zero, 
the smash sum A(BB has a deterministic scaling limit in W^. Before stating 
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our main results, we describe some related models and discuss our technique 
for identifying their common scaling limit, which comes from the theory of 
free boundary problems in PDE. 

The Diaconis-Fulton smash sum generalizes the model of internal diffu- 
sion limited aggregation ("internal DLA") studied in [20], and in fact was 
part of the original motivation for that paper. In classical internal DLA, 
we start with n particles at the origin o ^ 'L'^ and let each perform simple 
random walk until it reaches an unoccupied site. The resulting random set 
of n occupied sites in IJ^ can be described as the n-fold smash sum of {o} 
with itself. We will use the term internal DLA to refer to the same process 
run from an arbitrary starting configuration of particles. In this broader 
sense of the term, both the Diaconis-Fulton sum and the model studied in 
|20j are particular cases of internal DLA. 

In defining the smash sum A® various alternatives to random walk 
are possible. Rotor-router walk is a deterministic analogue of random walk, 
first studied by Priezzhev et al. j23] under the name "Eulerian walkers." At 
each site in 7? is a rotor pointing north, south, east or west. A particle 
performs a nearest-neighbor walk on the lattice according to the following 
rule: during each time step, the rotor at the particle's current location is 
rotated clockwise by 90 degrees, and the particle takes a step in the direction 
of the newly rotated rotor. In higher dimensions, the model can be defined 
analogously by repeatedly cycling the rotors through an ordering of the 2d 
cardinal directions in TJ^. The sum of two squares in 1? using rotor-router 
walk is pictured in Figure [Tj all rotors began pointing west. The shading 
in the figure indicates the final rotor directions, with four different shades 
corresponding to the four possible directions. 

The divisible sandpile model uses continuous amounts of mass in place 
of discrete particles. A lattice site is full if it has mass at least 1. Any 
full site can topple by keeping mass 1 for itself and distributing the excess 
mass equally among its neighbors. At each time step, we choose a full site 
and topple it. As time goes to infinity, provided each full site is eventually 
toppled, the mass approaches a limiting distribution in which each site has 
mass < 1; this is proved in [21]. Note that individual topplings do not 
commute. However, the divisible sandpile is "abelian" in the sense that any 
sequence of topplings produces the same limiting mass distribution; this is 
proved in Lemma |3.1[ Figure [T] shows the limiting domain of occupied sites 
resulting from starting mass 1 on each of two squares in Z^, and mass 2 on 
the smaller square where they intersect. 

Figure [T] raises a few natural questions: as the underlying lattice spacing 
becomes finer and finer, will the smash sum A (B B tend to some limiting 
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shape in M'^, and if so, what is this shape? Will it be the same limiting 
shape for all three models? To see how we might identify the limiting shape, 
consider the divisible sandpile odometer function 

u{x) = total mass emitted from x. 

Since each neighbor y ^ x emits an equal amount of mass to each of its 2d 
neighbors, the total mass received by x from its neighbors is 2^ '^Zy^x ^(?/)) 
hence 

/\u{x) = v{x) - a{x) (1) 

where (j{x) and ^{x) are the initial and final amounts of mass at x, re- 
spectively. Here A is the discrete Laplacian in Z"^, defined by An(x) = 

Equation ([T]) suggests the following approach to finding the limiting 
shape. We first construct a function on whose Laplacian is a — 1; an 
example is the function 

7(x) = - ^ 5ri(x,y)c7(y) (2) 

where in dimension d > 2> the Green's function gi{x,y) is the expected 
number of times a simple random walk started at x visits y (in dimension 
d = 2 we use the recurrent potential kernel in place of the Green's function). 
The sum u+j is then a superharmonic function on Z*^; that is, A{u+^) < 0. 
Moreover if / > 7 is any superharmonic function lying above 7, then f—j—u 
is superharmonic on the domain D = {x £ U^\y{x) = 1} of fully occupied 
sites, and nonnegative outside -D, hence nonnegative everywhere. Thus we 
have proved the following lemma of [21J. 

Lemma 1.1. Let a he a nonnegative junction on Z'^ with finite support. 
Then the odometer function for the divisible sandpile started with mass a{x) 
at each site x is given by 

M = S — 7 

where 7 is given by and 

s(x) = inf{/(x)|/ is superharmonic on Z"^ and / > 7} 
is the least superharmonic majorant of 7. 

For a reformulation of this lemma as a "least action principle," see 
Lemma 13.21 
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Figure 2: The obstacles 7 corresponding to starting mass 1 on each of two 
overlapping disks (top) and mass 100 on each of two nonoverlapping disks. 

Lemma [l . 1 1 allows us to formulate the problem in a way which translates 
naturally to the continuum. Given a function a on M.'^ representing the 
initial mass density, by analogy with Q we define the obstacle 

7(2;) = -|x|^ - / g{x,y)a{y)dy 

where g{x, y) is the harmonic potential on proportional to \x — yp"'^ in 
dimensions d > 3 and to — log |x — y| in dimension two. We then let 

s{x) = inf{/(x)|/ is continuous, superharmonic and / > 7}. 

The odometer function for a is then given by u = s — 7, and the final domain 
of occupied sites is given by 

D = {xe W^\s{x) > j{x)}. (3) 

This domain D is called the noncoincidence set for the obstacle problem with 
obstacle 7; for an in-depth discussion of the obstacle problem, see pO] • 

If A, B are bounded open sets in M'^, we define the smash sum of A and 
B as 

A®B = AlJBlJD (4) 
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where D is given by ([s]) with a = 1a + ^b- In the two-dimensional setting, 
an alternative definition of the smash sum in terms of quadrature identities 
is mentioned in |14] . 

By analogy with the discrete case, we would expect that 

C{A®B)=C{A) + C{B), 



where £, denotes Lebesgue measure in M"^. This is proved in Corollary 2.13 
More generally, if is a superharmonic function on Z"^, and ct is a mass 
configuration for the divisible sandpile (so each site j; G Z*^ has mass 
the sum ^^^.g^d h{x)a{x) can only decrease when we perform a toppling. 
Thus 

^ h{x)u{x) < h{x)a{x), (5) 

where is the final mass configuration. We therefore expect the domain D 
given by ^ to satisfy the quadrature inequality 



h{x)dx < I h{x)a{x)dx (6) 
D Jd 

for all integrable superharmonic functions h on D. For a proof under suitable 



smoothness assumptions on a and h, see Proposition 2.11 see also [27]. 

A domain Z) C M'^ satisfying an inequality of the form ^ is called a 
quadrature domain for a. Such domains are widely studied in potential 
theory and have a variety of applications in fluid dynamics [6l [25] . For more 
on quadrature domains and their connection with the obstacle problem, see 
[UEIIIHIIISIEZIISH]- Equation ([5| can be regarded as a discrete analogue of a 
quadrature inequality; in this sense, the three aggregation models studied in 
this paper produce discrete analogues of quadrature domains. Indeed, these 
aggregation models can be interpreted as discrete analogues of Hele-Shaw 
flow |6l [31] , which produces quadrature domains in the continuum. 

The main goal of this paper is to prove that if any of our three aggrega- 
tion models - internal DLA, rotor-router, or divisible sandpile - is run on 
finer and finer lattices with initial mass densities converging in an appropri- 
ate sense to a, the resulting domains of occupied sites will converge in an 
appropriate sense to the domain D given by Q. 

In order to state our main result, let us define the appropriate notion 
of convergence of domains, which amounts essentially to convergence in the 
Hausdorff metric. Fix a sequence 5n j representing the lattice spacing. 
Given domains An C (^nZ°' and D C M'^, write An — > D if for any e > 

D, n JnZ'^ CAnCD' (7) 
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for all sufficiently large n. Here and D'^ denote respectively the inner 
and outer e-neighborhoods of D: 



De = {x£ D\B{x,e) C D} 



..... (8) 



where B{x,€) is the ball of radius e centered at x. 

For X G JnZ'^ we write x^ = [x+ ^,x - ^] . For t G M write [t] for 
the closest integer to t, rounding up if t G Z + ^. 

The initial data for each of our lattice models consists of a function 
an '■ ^rJ^^ — > Z>o representing the number of particles (or amount of mass) 
at each lattice site; we think of density with respect to counting 

measure on (J^Z'^, and refer to it as the "initial density." 

Throughout this paper, to avoid trivialities we work in dimension d >2. 
Our main result is the following. 

Theorem 1.2. Let n C R'^ , d > 2 be a bounded open set, and let a : M — > 
Z>o be a bounded function which is continuous almost everywhere, satisfying 
{(7 > 1} = 0,. Let Dn,Rn,In be the domains of occupied sites formed from 
the divisible sandpile, rotor-router model, and internal DLA, respectively, in 
the lattice 6n'^'^ started from initial density 



(T„(x) 

Then as n ] oo 



Dn,Rn Dun- 

and if 6n log n [ 0, then with probability one 

In ^ Dun 

where D is given by and the convergence is in the sense of 

Remark. When forming the rotor-router domains Rn, the initial rotors in 
each lattice 5„Z'^ may be chosen arbitrarily. 

As an immediate consequence, all three lattice models have a rotationally- 
invariant scaling limit. This follows from the rotational symmetry of the 
obstacle problem: if p is a rotation of M'^, and D is the noncoincidence set 
^ for initial density a, then the noncoincidence set for initial density aop~^ 
is pD. 
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We prove a somewhat more general form of Theorem 1.2 which allows 
for some flexibility in how the discrete density cr„ is constructed from a. 
In particular, taking cr = 1^ + Ig we obtain the following theorem, which 
explains the similarity of the three smash sums pictured in Figure [T| 

Theorem 1.3. Let A,BcM.'^ be bounded open sets whose boundaries have 
measure zero. Let Dn, Rn, In be the smash sum of An br^^ and B n 6n'^'^, 
formed using divisible sandpile, rotor-router and internal DLA dynamics, 
respectively. Then as n ] oo 

Dn,Rn -^A®B; 
and if 5n log | 0, then with probability one 

In^A®B 

where A® B is given by and the convergence is in the sense of 



For the divisible sandpile. Theorem 1.2 can be generalized by dropping 



the requirement that a be integer valued; see Theorem 3.9 for the precise 
statement. Taking a real- valued is more problematic in the case of the 
rotor-router model and internal DLA, since these models work with discrete 
particles. Still, one might wonder if, for example, given a domain A C 
W^, starting each even site in A n 5nZ°' with one particle and each odd 
site with two particles, the resulting domains Rn,In would converge to the 
noncoincidence set D for density a = |1a- This is in fact the case: if an is 
a density on Sn*^^ J as long as a certain "smoothing" of Cn converges to (j, 
the rotor-router and internal DLA domains started from initial density cr„ 



will converge to D. See Theorems 4.8 and 5.1 for the precise statements. 

One interesting case not covered by Theorems |1.2| and |1.3| is the case of 
multiple point sources. Lawler, Bramson and Griffeath [20 showed that the 
scaling limit of internal DLA in Z'^ with a single point source of particles is 
a Euclidean ball. In [21] , the present authors showed that for rotor- router 
aggregation and the divisible sandpile in Z*^ with a single point source, the 
scaling limit in both cases is also a Euclidean ball. 

For X G M'^ write x" for the closest lattice point in 6fi^'^ .} breaking ties to 
the right. Our shape theorem for multiple point sources, which is deduced 
from Theorem 1.3 using the main results of [20] and [21], is the following. 

Theorem 1.4. Fix xi, . . . ,Xk G M"^ and Ai, . . . , > 0. Let Bi be the ball 
of volume Aj centered at Xj. Fix a sequence 6n i 0, and for x G (5^^°' let 

k 



an{x) 



{x=xv} 
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Figure 3: The rotor-router model in 1? started from two point sources on 
the X-axis. The boundary of the limiting shape is an algebraic curve (11 ) of 
degree 4. 



Let Dn, Rn, In be the domains of occupied sites in 6n'^'^ formed from the di- 
visible sandpile, rotor-router model, and internal DLA, respectively, started 
from initial density Un- Then as n —> oo 

Dn,Rn^ Bi®...®Bk\ (9) 

and if 5n <^ /n, then with probability one 

In^ Bi® ...®Bk 

where © denotes the smash sum and the convergence is in the sense 

Implicit in ^ is the associativity of the smash sum operation, which is 
not readily apparent from the definition Q. For a proof of associativity, see 



Lemma 6.1, Alternatively, the associativity follows from basic properties of 
"partial balayage" operators on measures on M*^; see Theorem 2.2(iii)]. 
For related results in dimension two, see |261 Prop. 3.10] and [31, section 
2.4]. 

We remark that a similar model of internal DLA with multiple point 
sources was studied by Gravner and Quastel [12] , who also obtained a vari- 



ational solution. In their model, instead of starting with a fixed number of 
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particles, each source Xi emits particles according to a Poisson process. The 
shape theorems of [12] concern convergence in the sense of volume, which is 
a weaker form of convergence than that studied in the present paper. 

In Proposition 6.6 we show that the smash sum of balls i?i © . . . © Bf^ 
arising in Theorem 1.4 obeys the classical quadrature inequality 



k 

h{x)dx <'^Xih{xi) (10) 



for all integrable superharmonic functions h on Si © ... © i?^, with equality 
if h is harmonic. This can be regarded as a generalization of the classical 
mean value property of harmonic functions, which corresponds to the case 
A; = 1. 

Let Bi, . . . , Bj^ be disks in with distinct centers xi, . . . , x^. From m\ 
and work of Aharonov and Shapiro [T and Gustafsson |13l ITS] it follows that 
that the boundary d[Bi © . . . © i?fc) lies on an algebraic curve of degree 2k. 
More precisely, there is a polynomial P G M[x, y] of the form 

P(x, y) = (x^ + y^)^ + lower order terms 

and there is a finite set of points E possibly empty, such that 

d{Bi ®...(BBk) = {{x, y) G R^\P{x, y) = 0} - E. 

For example, if Bi and B2 are disks of equal radius r > 1 centered at 
(1,0) and (—1,0), then d{Bi © B2) is given by the quartic curve |29] 

(x2 + 2/2)2 _ (^2 ^ y2^ _ 2(^2 _ y2^ ^ Q (11) 

This curve describes the shape of the rotor-router model with two point 
sources pictured in Figure [3| 



2 Potential Theory Background 

In this section we review the basic properties of superharmonic potentials 
and of the least superharmonic majorant. The proofs are deferred to sec- 
tion [7j In view of the diverse readership who might be interested in the 
models we study, we have elected to include somewhat more than the usual 
amount of background material. 
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2.1 Least Superharmonic Majorant 

Since we will often be working with functions on M'^ which may not be twice 
differentiable, it is desirable to define superharmonicity without making ref- 
erence to the Laplacian. Instead we use the mean value property. A function 
u on an open set $7 C M'^ is superharmonic if it is lower-semicontinuous and 
for any ball B(x, r) C O 

u{x) > Aru{x) := I u{y)dy. (12) 

^dr"^ JB{x,r) 

Here uJd is the volume of the unit ball in M'^. We say that u is subharmonic 
if —u is superharmonic, and harmonic if it is both super- and subharmonic. 

The following properties of superharmonic functions are well known; for 
proofs, see e.g. |3], [8] or [22] . 

Lemma 2.1. Let u be a superharmonic function on an open set 17 C M*^ 
extending continuously to 0,. Then 

(i) u attains its minimum in 0, on the boundary. 

(a) Ifhis continuous on 0,, harmonic on 0,, and h = u on dO,, then u > h. 
(Hi) If B{x,ro) C B{x,ri) C Jl, then 

Arf^u{x) > Arj^u{x). 



(iv) If u is twice differentiable on 0,, then Au < on 0,. 

(v) If B C 0, is an open ball, and v is a function on Q which is har- 
monic on B, continuous on B, and agrees with u on B'^, then v is 
superharmonic. 

Given a function 7 on which is bounded above, the least superhar- 
monic majorant of 7 (also called the solution to the obstacle problem with 
obstacle 7) is the function 

s{x) = inf{/(x)|/ is continuous, superharmonic and / > 7}. (13) 

Note that since 7 is bounded above, the infimum is taken over a nonempty 
set. 

Lemma 2.2. Let -f be a uniformly continuous function which is bounded 
above, and let s be given by (13). Then 
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(i) s is superharmonic. 

(ii) s is continuous. 

(Hi) s is harmonic on the domain 

D = {xe M^|s(x) > -fix)}. 



2.2 Superharmonic Potentials 

Next we describe the particular class of obstacles which relate to the aggre- 
gation models we are studying. For a bounded measurable function a on M.'^ 
with compact support, write 



Ga{x) 



/ g{x,y)a{y)dy, (14) 



where 



-|log|x-y|, d = 2; 
ad\x-y\^-'^, d>3. 



9ix,y) = { : as) 



Here = where ojd is the volume of the unit ball in M . Note that 

(15) differs by a factor of 2d from the standard harmonic potential in W^; 
however, the normalization we have chosen is most convenient when working 
with the discrete Laplacian and random walk. 

The following result is standard; see |8j Theorem 1.1.7.2]. 

Lemma 2.3. Let a he a measurable function on with compact support. 

(i) If a is bounded, then Ga is continuous. 

(ii) If a is C^, then Ga is G^ and 

AGa = -2da. (16) 



Regarding (ii), if we remove the smoothness assumption on a, equation 
(161 remains true in the sense of distributions. For our applications, however. 



we will not need this, and the following lemma will suffice. 

Lemma 2.4. Let a be a bounded measurable function on Mf^ with compact 
support. If a > on an open set 17 C M"', then Ga is superharmonic on 17. 
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By applying Lemma 2.4 both to a and to —a, we obtain that Ga is harmonic 
off the support of a. 

Let B = B{o, r) be the ball of radius r centered at the origin in W^. We 
compute in dimensions d > 3 



Likewise in two dimensions 



2r2 / r V 
J-2\\x\) 



X > r. 



jr2(l - 21ogr) - |x| < r 

I — 2r log |x| > r. 

Fix a bounded nonnegative function a on with compact support, and 

let 

7(x) = - Gct(x). (19) 

Let 

s(x) = inf{/(x)|/ is continuous, superharmonic and / > 7} (20) 
be the least superharmonic majorant of 7, and let 

D = {x(^ W^\s{x) > 7(x)} (21) 
be the noncoincidence set. 

Lemma 2.5. (i) j{x) + |xp is subharmonic on M*^. 

(ii) If cr < M on an open set 17 C MJ^, then 7(x) — (M — l)|xp is super- 
harmonic on 0,. 



Lemma 2.6. Let u = s — j, where 7 and s are given by (19) and (20). 
Then 

(i) u{x) — |xp is superharmonic on W^. 

(ii) If a < M on an open set 17 C M*^, then u{x) + M|xp is subharmonic 
on 0,. 



Lemma 2.7. Let D be given by ^2^. Then {a > 1}° C D. 

The next lemma concerns the monotonicity of our model: by starting 
with more mass, we obtain a larger odometer and a larger noncoincidence 
set; see also [271 Lemma 7.1]. 
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Lemma 2.8. Let ai < cj2 be functions on M*^ with compact support, and 
suppose that a2{x)dx < oo. Let 

7i(x) = - GfTj(x), i = l,2. 

Let Si be the least superharmonic majorant of 'ji, let Ui = Si — 'ji, and let 

Di = {x\si{x) > 7j(x)}. 

Then ui < U2 and Di C D2. 

Our next lemma shows that we can restrict to a domain 17 C M*^ when 
taking the least superharmonic majorant, provided that 0, contains the non- 
coincidence set. 



Lemma 2.9. Let j,s,D be given by (19)-(21). Let Q. C be an open set 
with D C 0,. Then 

s{x) = inf{/(a;)|/ is superharmonic on continuous, and f > 7}. (22) 
2.3 Boundary Regularity for the Obstacle Problem 



Next we turn to the regularity of the solution to the obstacle problem (20) 
and of the free boundary dD. There is a substantial literature on boundary 
regularity for the obstacle problem. In our setting, however, extra care must 



be taken near points where (t{x) = 1: at these points the obstacle (19) is 
harmonic, and the free boundary can be badly behaved. 

One approach, adopted in [17J, is to reformulate the obstacle problem 
in M'^ as an obstacle problem on the set {a > 1}. The results of [17j apply in 
the situation that a > 1 everywhere, but we will want to relax this condition. 
We show only the minimal amount of regularity required for the proofs of 
our main theorems. Much stronger regularity results are known in related 
settings; see, for example, jHIS]. 

The following lemma shows that if the obstacle is sufficiently smooth, 
then the superharmonic majorant cannot separate too quickly from the ob- 
stacle near the boundary of the noncoincidence set. As usual, we write 

D^ = {x e D\B{x,e) C D} (23) 

for the inner e-neighborhood of D. 

Lemma 2.10. Let a be a function on M'^ with compact support. Let 



7,8,1? be given by (19)-(21), and write u = s — ■j. Then u is C , and for 



y G dDe we have \Vu{y)\ < Coe, for a constant Co depending on a. 
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Next we show that mass is conserved in our modeh the amount of mass 
starting in D is Jj^ a{x)dx, while the amount of mass ending in D is C{D), 
the Lebesgue measure of D. Since no mass moves outside of D, we expect 
these two quantities to be equal. Although this seems intuitively obvious, the 
proof takes some work because we have no a priori control of the domain D. 
In particular, we first need to show that the boundary of D cannot have 
positive d- dimensional Lebesgue measure. 

Proposition 2.11. Let a be a function on with compact support, 
such that £(cr^^(l)) = 0. Let D he given by (21). Then 



(i) C{dD) = 0. 

(ii) For any function h G C^{D) which is superharmonic on D, 



h{x)dx < / h{x)a{x)dx. 
D Jd 

Note that by applying (ii) both to h and to —h, equality holds whenever 
h is harmonic on D. In particular, taking h = 1 yields the conservation of 
mass: J"^ a{x)dx = C{D). 



We will also need a version of Proposition 2.11 which does not assume 
that a is or even continuous. We can replace the assumption and the 
condition that £(cr^^(l)) = by the following condition. 

For ah X gR'^ either a{x) < A or a{x) > 1 (24) 

for a constant A < 1. Then we have the following result. 

Proposition 2.12. Let a be a bounded function onM'^ with compact support. 



Let D be given by (21), and let D = Dujo" > 1}°. If a is continuous almost 



everywhere and satisfies (24), then 

(i) C{db) = 0. 

(ii) >C(L') = f^a{x)dx. 

(Hi) C[D) = jj^cr{x)dx. 

In particular, taking cr = 1^ + 1b, we have 

A®B = AiJB\JD = b 

where ® denotes the smash sum (|4]). From (iii) we obtain the volume addi- 
tivity of the smash sum. 
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Corollary 2.13. Let A,BcM. be bounded open sets whose boundaries 
have measure zero. Then 

C{A®B) =C{A) + C{B). 

Our next lemma describes the domain resulting from an initial mass 
density m > 1 on a ball in W^. Not surprisingly, the result is another ball, 
concentric with the original, of m times the volume. In particular, if m is 
an integer, the m-fold smash sum of a ball with itself is again a ball. 



Lemma 2.14. Fix m > 1, and let D be given by (21) with a = mlB[o,r)- 
Then D = B{o,m^/'^r). 

Proof. Since 7(x) = — |xp — Ga{x) is spherically symmetric and the least 
superharmonic majorant commutes with rotations, D is a ball centered at 
the origin. By Proposition 2.12[ ii) we have C{D) = m£{B{o,r)). □ 



Next we show that the noncoincidence set is bounded; see also |27| 
Lemma 7.2]. 

Lemma 2.15. Let a be a function on with compact support, satisfying 



< a < M. Let D be given by (21). Then D is bounded. 



Pro of. Le t B = B{o, r) be a ball containing the support of a. By Lemmas 2.8 
and 



2.14 



the difference s — 7 is supported in B{o, M^^^r). □ 



2.4 Convergence of Obstacles, Majorants and Domains 

In broad outline, many of our arguments have the following basic structure: 

convergence of densities =^ convergence of obstacles 

=^ convergence of majorants 
=^ convergence of domains. 

An appropriate convergence of starting densities is built into the hypotheses 
of the theorems. From these hypotheses we use Green's function estimates 
to deduce the relevant convergence of obstacles. Next, properties of the 
obstacle can often be parlayed into corresponding properties of the least 
superharmonic majorant. Finally, deducing convergence of domains (i.e., 
noncoincidence sets) from the convergence of the majorants often requires 
additional regularity assumptions. The following lemma illustrates this basic 
three-step approach. 
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Lemma 2.16. Let a and an, n = 1,2, . . . be densities on M'^ satisfying 

< cr,cr„ < MIb 
for a constant M and ball B = B{o, R). If an ^ cr in , then 
(i) Gan — > Ga uniformly on compact subsets ofW^. 

(a) s„ — > s uniformly on compact subsets o/M"', where s,s„ are the least 
superharmonic majorants of the functions 7 = —\x\'^ — Ga, and 7^ = 
— |xp — Gan, respectively. 

(Hi) For any e > we have C D(^n) for all sufficiently large n, where 
D,D^n) o,re the noncoincidence sets {s > 7} and {s„ > 7„}, respec- 
tively. 

According to the following lemma, a nonnegative function on with 
bounded Laplacian can grow at most quadratically. 

Lemma 2.17. Fix < (3 < 1. There is a constant cp such that any non- 
negative function f on with f{o) = and |A/| < X in B{o,R) satisfies 

/(x) < c/3A|x|^ x£B{o,pR). 



We will also use the following continuous version of Lemma |2.17 



Lemma 2.18. Fix < (3 < 1. There is a constant c'^ such that the following 
holds. Let f be a nonnegative function on R"' with f{o) = 0, and let 



f±{x) = f{x)±\^^. 

If f+ is subharmonic and /_ is superharmonic in B(o,R), then 

/(x) < c'^A|x|2, xeB{o,f3R). 

In the following lemma, denotes the inner e-neighborhood of A, as 
defined by ([23|. 

Lemma 2.19. Let AjBcW^ be bounded open sets. For any e > there 
exists rj > with 

{AUB)^ C ArjUBr,. 
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points 


x' = 








sets 


= 


A n (5„Z<^ 




functions 


r = 






= fix----) 



Table 1: Notation for transitioning between Euclidean space and the lattice. 



2.5 Discrete Potential Theory 

Fix a sequence Sn [ with Si = 1. In this section we relate discrete su- 
perharmonic potentials in the lattice Sn'^'^ to their continuous counterparts 
in W^. If A is a domain in ^rZ*^, write 



A° = A- 



for the corresponding union of cubes in W^. If A is a domain in W^, write 
A- = Ar] 5n1''^. Given xeR'^, write 



— ,x+ - 



for the closest lattice point to x, breaking ties to the right. For a function 
/ on Sn'^'^, write 

f°{x) = fix----) 

for the corresponding step function on M''. Likewise, for a function / on M"^, 
write /•■ = /i^^^d. These notations are summarized in Table 1. 
We define the discrete Laplacian of a function / on Snl''^ to be 



According to the following lemma, if / is defined on M and is sufficiently 
smooth, then its discrete Laplacian on c^^Z'^ approximates its Laplacian 
on W^. 

Lemma 2.20. // / : M*^ — ^ M has continuous third derivative in a 6n- 
neighborhood of x £ , and A is the maximum 'pure third partial of f in 
this neighborhood, then 
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In three and higher dimensions, for x, y G b^JH^ write 

(7„(x, y) = bl~^gi{b~^x, d'^), d>3, (25) 

where 

g,{x,y)=E^#{k\Xk = y} 



is the Green's function for simple random walk on Z . The scaling in (25) 
is chosen so that Axgn{x,y) = —6~'^l^x=yy In two dimensions, write 

9n{x, y) = -a{5^ x,6n y) + -log5n, d = 2, (26) 

TT 

where 

a{x, y) = lim ( E^; < m\Xk = x} - E^.#{A; < ■m\Xk = y}) 

m— >oo 

is the recurrent potential kernel on 1?. 
Lemma 2.21. In all dimensions d >2, 

9n{x, y) = g{x, y) + 0{5l\x - y\-'^) 



where g is given by (15). 

Our next result adapts Lemma 2.16[ i) to the discrete setting. We list 
here our standing assumptions on the starting densities. Let a be a function 
on with compact support, such that 

< 0- < M (27) 

for some absolute constant M. Suppose that a satisfies 

C{DC{a))=G (28) 

where DC {a) denotes the set of points in where a is discontinuous. 
For n = 1, 2, . . . let (Tn be a function on Sn'^'^ satisfying 

< o-„ < M (29) 

and 

(7°(x) ^ cj(x), x i DC{a). (30) 

Finally, suppose that 



There is a ball i? C M*^ containing the supports of a and for all n. 



(31) 
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Although for maximum generality we have chosen to state our hypotheses 
on a and an separately, we remark that the above hypotheses on (T„ are 
satisfied in the particular case when (t„ is given by averaging a over a small 
box: 



an{x) = 



a{y)dy. 



In parallel to (14), for x S 5^^"^ write 



Gn(yn{x) = 5i ^ gn{x,y)an{y), 



(32) 



(33) 



where gn is given by (25) and (26). Since gn{x, •) is discrete harmonic except 
at X, where it has discrete Laplacian we have for any function / 

on with bounded support 



(34) 



The first equality can be seen directly from the symmetry gn{x, y) = gniu, x) 
of the Green's function, or else by observing that the operators A and Gn 
are both defined by convolutions, so they commute. 



Lemma 2.22. If a, an satisfy (21)-(31), then 

{Gncrnf^ Ga 
uniformly on compact subsets ofW^. 



Our next result adapts Lemma 2.9 to the discrete setting. Let an be a 
function on 6n'^'^ with finite support, and let 7„ be the function on 
defined by 



ln{x) 



\X\ 



Gnan{x). 



(35) 



Let 



Sn = inf{/(2;)|/ is superharmonic on (5nZ'^ and / > 7„} (36) 
be the discrete least superharmonic majorant of 7„, and let 



Dn = {xe 6n Z \Sn{x) > Jn{x)} . 



(37) 



Lemma 2.23. Let jn,Sn,Dn be given by {35)-{37). If ^ C Sn'^'^ satisfies 
Dn C then 



Sn{x) = inf{/(x)|/ is superharmonic on 0, and / > 7^}. 



(38) 
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3 Divisible Sandpile 



3.1 Abelian Property and Least Action Principle 

It will be useful here to prove the abelian property for the divisible sandpile 
in somewhat greater generality than that of [21j. We work in continuous 
time. Fix r > 0, and let T : [0, r] — > Z'^ be a function having only finitely 
many discontinuities in the interval [0, t] for every t < t. The value T(t) 
represents the site being toppled at time t. The odometer function at time 
t is given by 

ut{x)=C{T-\x)n[0,t]) , 

where C denotes Lebesgue measure. We will say that T is a legal toppling 
function for an initial configuration r]Q if for every < t < r 

vtim) > 1 

where 

Vtix) = r]o{x) + Aut{x) 

is the amount of mass present at x at time t. If in addition rjr < 1, we say 
that T is complete. 

Lemma 3.1. (Abelian Property) If Ti : [0,ti] and T2 : [0,T2\ IJ^ 

are complete legal toppling functions for an initial configuration r]Q, then for 
any x G Z*^ 

C {T,-Hx)) = C {T,\x)) . 
In particular, ti = T2 and the final configurations rj^^ , r/^-j are identical. 
Proof. For i = 1,2 write 

ui{x)=C{Tr\x)n[0,t]); 
4{x) = rio{x)+Aui{x). 

Write u' = and rj' = r]\^. Let t(0) = and let t(l) < i(2) < ... be 
the points of discontinuity for Ti. Let x^ be the value of Ti on the interval 
{t{k — l),t{k)). We will show by induction on k that 

u^{xk) > ^^t(fc)(a;fe). (39) 

Note that for any a; 7^ x^, if u^(^i^^{x) > 0, then letting j < A; be maximal 
such that Xj = x, since Ti 7^ x on the interval (t{j),t(k)), it follows from 
the inductive hypothesis ( 39 ) that 

= u'^{xj) > = ul^^^{x). (40) 
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Since Ti is legal and T2 is complete, we have 
hence 

Since Ti is constant on the interval (t{k — l),t(k)) we obtain 



By (40), each term in the sum on the right side is nonnegative, completing 
the inductive step. 



Since t{k) j ti as t 00, the right side of (40) converges to u^{x) as 
k I cxD, hence > u^. After interchanging the roles of Ti and T2, the result 
follows. □ 

The next lemma is a simple reformulation of Lemma but is some- 
times more convenient to use. 

Lemma 3.2. (Least Action Principle) Let a be a nonnegative function on Z*^ 
with finite support, and let u he its divisible sandpile odometer function. Lf 
ui : X'^ —>■ M>o satisfies 

cr + Ani<l (41) 

then ui > u. 

Proof. From Lemma |l.l| we have u = s — 7. Since 

Aui < 1 — o" = — A7 

the sum ui + 7 is superharmonic on U^. Since ui is nonnegative, it follows 
that Ml + 7 > s, hence ui > u. □ 

Note that ui need not have finite support. The moniker "least action 
principle" comes from the following interpretation. Starting from a, and 
toppling each site x G Z"^ for time ui{x) without regard to whether these 
topplings are legal, equation (41) says that the resulting configuration is 
stable: no site has mass greater than 1. Clearly the odometer u is such a 
stabilizing function. The Least Action Principle says that u is (pointwise) 
minimal among all such functions. 
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3.2 Convergence of Odometers 

By the odometer function for the divisible sandpile on br^'^ with initial 
density we will mean the function 

"UnC^;) = (^^-total mass emitted from x if each site y starts with mass (Tn{y)- 

Theorem 3.3. Let Un be the odometer function for the divisible sandpile 
on 6n'^'^ with initial density an- If cf, an satisfy ( 27)-( 31), then 

—>■ s — ^ uniformly, 

where 

■~f{x) = — |xp — Ga{x), 



Ga is given by (I4), and s is the least superharmonic majorant ofj. 

Lemma 3.4. Let Dn be the set of fully occupied sites for the divisible sand- 
pile in br{L'^ started from initial density an ■ There is a ball ncR'^ with 

\jD°UDcn. 

n>l 

Proof. Let An be the set of fully occupied sites for the divisible sandpile in 
JnZ'^ started from initial density t{x) = MIx^b, where B is given by (31). 
From the abelian property, Lemma |3.1[ we have D„ C An- By the inner 
bound of [21' Theorem 1.3] if we start with mass m = 25~'^C{B) at the 
origin in 5n'^'^, the resulting set of fully occupied sites contains B '-, by the 
abelian property it follows that if we start with mass Mm at the origin in 
5„Z'^, the resulting set i7„ of fully occupied sites contains An- By the outer 
bound of |2H Theorem 1.3], 0„ is contained in a ball O of volume 2>MC{B). 
By Lemma 2.15 we can enlarge O if necessary to contain D. □ 

For X G 5n'^'^ write 

7n(x) = - G„cr„(x), 

where On is defined by ( |33| ). Denote by Sn the least superharmonic majorant 
of 7„ in the lattice 6n'^ ■ 



Lemma 3.5. Let VL be as in Lemma 3.4 ■ There is a constant M' independent 
of n, such that |7„| < M' in Jl--. 
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Proof. By Lemma 2.21 we have for x G 



\Gn(Tn{x)\ <5^Y^ \gn{x,y)\an{y) 

<2M5t \9ix^y)\ 
< CMR^ log R 

for a constant C depending only on d, where R is the radius of O. It follows 
that |7„| < (CM + l)i?2logi2in J^^ □ 

Lemma 3.6. Fix x G M"^, and for y G (^„Z°' let 



a{y) = Sngnix'-'',y) - I g{x,z)dz. 



There are constants Ci,C2 depending only on d, such that 

(i) \a{y)\<Ci5l+''\x-y\^-''. 

(ii) Ifyi ~ 2/2, then \a{yi) - a[y2)\ < C25l'^'^\x - yi\-'^. 



(42) 



Proof, (i) By Lemma 2.21, if z G then 



\gn{x-',y) - g{x, z)\ < \gn{x-'-,y) - gix-'-,y)\ + y) - g{x, z)\ 

< C6n\x - y]^-"" 

for a constant C depending only on d. Integrating over y^ gives the result. 



(ii) Let p = y2 — yi- By Lemma 2.21 , we have 



a(yi) - a(y2)| < / \gn{x--,yi) - gnix--,y2) - gix,z) + g{x,z +p)\dz 



< / \9{x'',yi) - g{x-'-,y2) - g{x,z) + g{x,z +p)\dz + 
Jy? 

+ 0(<5^2|^-yi|-'^). 
Writing q = x'' — x — yi + z, the integrand can be expressed as 
\f{o)-f{p)-fiq)+f{p + q)\, 

where 

f{w) = g{w,x-- - yi). 
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Since \p\, \q\, \p + q\ < 5n(v^ + l), we have by Taylor's theorem with remain- 
der 



\Ko) - f{p) - m +f{p + q)\< M5l 



dxidxj 



< 



\24/TT)6l\x-yi\-^, d = 2 

6d^{d-l){d-2)ad6l\x-yi\-'^, d>3. □ 



Lemma 3.7. Let 17 be an open ball as in Lemma 3.4 ^ and let J7i be a ball 
with C fii- Let 

(j)n = -A(s„1q^:). 

Then 





uniformly on 0,. 



Proof. Let Un{x) be the amount of mass present at x in the final state of 
the divisible sandpile on S'li"^^ . By LeiniTia. |l - ll we have Sn — Un -\- 'yn-j tience 

As„ = {l^n - CTn) + (0"„ - 1) = - 1. 



In particular, |As„| < 1. 
Note that 

whereas Sn = Gn(t)n in O'', by (|34]). Hence for x G J7 

G(t)^{x) - s°(x) = - G„(/)„(x--) 

= - ^ (l)n{y)a{y) 



(43) 



where a{y) is given by (42). Now write 

4)n = -{/\Sn)ln--l + /l - /2 

where /i is supported on the inner boundary of VL'l and given there by 
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while /2 is supported on the outer boundary of and given there by 

1 



f2{z) 



The sum in (|43|) now takes the form 



As„(y)a(y) + ^ 



Sn{y)a{z) - Sn{z)a{y) 



y^z 



2d5l 



(44) 

Let i?, ^1 be the radii of By Lemma |3.6[ i), summing in spherical 

unded ir 

{Ci5i+'^r^-'^)dr = CiduJdRiS„ 



shells about x, the first sum in (44) is bounded in absolute value by 



To bound the second sum in (44), note that s„ = 7n outside so 

\sn{y)oi{z) - Sn{z)a{y)\ < \jn{y)\\a{y) -a{z)\ + \a{y)\\'jn{y) -ln{z)\. (45) 

By Lemmas |3.5| and |3.6[ ii), the first term is bounded by 

hn{y)\Hy) - «(^)l < C2M'6l"^''\x - y\-''. 

Fix e > 0, and let O2 be a ball with Cli C O2. Since 7 is uniformly continuous 
on O2, and 7n — > 7 uniformly on O2 by Lemma 2.22[ for sufficiently large n 
we have 

hn{y) - ln{z)\ < hn{y) - 7(^)1 + h{y) - l{z) \ + h{z) - -YnizM < €. 



Thus by Lemma |3.6[ i) the second term in (45) is bounded by 

Hy)\\ln{y) - ln{z)\ < Ci6i+''\x - 



e. 

Since x € 0, and y is adjacent to d^li, we have |x — y| > Ri — R — 6n, so the 
second sum in (44) is bounded in absolute value by 

6~^#dn'i {c2M'6l+\Ri - R)-^ + Ci5l+'^{Ri - R)^-'^e) < C^e 

for sufficiently large n. □ 
Lemma 3.8. — > s uniformly on compact subsets ofW^. 
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Proof. By Lemma 3.4 there is a ball 0, containing D and for all n. 
Outside 0, we have 



„□ 



,□ 



7 



uniformly on compact sets by Lemma 2.22 
To show convergence in Q, write 



s{x) 



x-y 
A 



(46) 



where 77 is the standard smooth mollifier 
rj{x) - 



Cei/(N^-i), \x\ < 1 
0, \x\ > 1 



normalized so that J^^^ rjdx = 1 (see |9l Appendix C.4]). Then s is smooth 
and superharmonic. Fix e > 0. By Lemma |2.2[ ii) and compactness, s is 
uniformly continuous on 0,, so taking A sufficiently small in (46) we have 
|s — s| < e in 17. Let be the maximum third partial of s in 0,. By 
Lemma I2. 201 the function 



qn{x) = s-'-{x) 



A.^6n\x\ 



is superharmonic in 0, '. By Lemma 2.22 we have 7^ 7 uniformly in 0, 



Taking N large enough so that gA^fJ^Jxp < e in 17 and |7„ — 7''| < e in 
for all n > A^, we obtain 

Qn > s" — e > s" — 2e > 7" — 2e > 7„ — 3e 

in In particular, the function = max(g„ + 3e, 7„) is superharmonic 
in 0, '. By Lemma 2.23 it follows that fn > Sn, hence 

■Sn < 9ra + 3e < s" + 3e < s" + 4e 

in Q". By the uniform continuity of s on 0, taking A larger if necessary we 
have Is — s''^| < e in $7, and hence sE^ < s + 5e in i7 for all n > N. 



For the reverse inequality, let 



where is an open ball containing J7. By Lemma 3.7 we have 



□ 



Gct>\ 



□ 



< e 
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and hence 



G0°>7:^-e>7-2e. 



(47) 



on 17 for sufficiently large n. Since (f>^ is nonnegative on fl, by Lemma 
the function G(j)^ is superharmonic on 17, so by (47) the function 



2.4 



ipn = max(G(?:)° + 2e, 7) 



is superharmonic on O for sufficiently large n. By Lemma 2.9 it follows that 
■^n > hence 

s^>G<P°-e>s-3e 
on n for sufficiently large n. □ 



1.1 



— > 7 and ^ s uniformly on $7. By Lemma 
Un = Sn — ^n- Since = = s — 7 off we conclude that u}^ 
uniformly. 

3.3 Convergence of Domains 



Proof of Theorem \3.S\ Let J7 be as in Lemma |3.4[ By Lemma s |2 . 22 1 and |3 .8 
we have 7^ 



we have 
^ s — 7 

□ 



In addition to the conditions (27)-(31) assumed in the previous section, we 

(48) 



assume in this section that the initial density a satisfies 

For ah X G M'^ either a{x) < A or ct(x) > 1 
for a constant A < 1. We also assume that 



{cj > 1} = {cr > 1}". (49) 
Moreover, we assume that for any e > there exists A'^(e) such that 

If X e {ct > 1},, then (t„(x) > 1 for all n > A^(e); (50) 

and 



If X ^ {cr > 1}', then (T„(x) < A for all n > N{e) 



(51) 



As before, we have chosen to state the hypotheses on a and o"„ separately for 
maximum generality, but all hypotheses on it„ are satisfied in the particular 



case when cr„ is given by averaging a in a small box ( 32 ) . 
We set 

7(x) = — |x|^ — Ga{x) 
with s the least superharmonic majorant of 7 and 

D = {xe M.'^\s{x) > 7(x)}. 
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We also write 

D = DU{xe R'^\a{x) > 1}°. 

For a domain yl C M'^, denote by and A"" its inner and outer open e- 
neighborhoods, respectively. 



Theorem 3.9. Let a and an satisfy (21)-(31) and (48)-(51). For n > 1 let 
Dn he the domain of fully occupied sites for the divisible sandpile in br^^ 
started from initial density an- For any e > we have for large enough n 



(52) 



According to the following lemma, near any occupied site x G lying 
outside -D*^, we can find a site y where the odometer Un is relatively large. 
This is a discrete version of a standard argument for the obstacle problem; 
see for example Friedman flO], Ch. 2 Lemma 3.1. 

Lemma 3.10. Fix e > and x € Z)„ with x ^ D^. If n is sufficiently large, 
there is a point y E Jn^*^ with \x — y\ < | + 5n O'lT-d 

Un{y) > Un{x) + ^—^e'^. 



Proof By (|49|, we have {fi >l} C D. Thus if x ^ D', the bah 5 = B{x, ^y- 
is disjoint from {a > 1}'^/^. In particular, if n > iV(§), then by (51 ) we have 
o"n < A on Thus the function 

w{y) = Un{y) - {I - \)\x- 

is subharmonic on BCiDn, so it attains its maximum on the boundary. Since 
w{x) > 0, the maximum cannot be attained on dDn, where Un vanishes; so 
it is attained at some point y G dB, and 

un{y) > w{y) + {1 - X) (^^y . 

Since w{y) > 'w{x) = Un{x), the proof is complete. □ 



Proof of Theorem \3.S\ Fix e > 0. By Lemma 2.19 we have 

b^C Dr^[j{a> l]r, 
for some ?? > 0. Let u = s — 7. Since the closure of is compact and 



contained in D, we have u > iriri on Drj for some m-q > 0. By Theorem 3.3 
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for sufficiently large n we have Un> u" 



i^m^ > on 



hence D]l C Dn- 



Likewise, by (50 ) we have {a > 1}^ C Dn for large enough n. Thus D^ C Dn- 
For the other inclusion, fix x G_^Z°' with x ^ D'^. Since u vanishes on 
the ball B = -B(x, §), by Theorem 



3.3 



we have Un < 



l-A^2 



on B--UdB-- for 

all sufficiently large n. By Lemma 3.10 we conclude that x ^ Dn, and hence 
Dn C D'--- □ 



4 Rotor-Router Model 



In trying to adapt the proofs of Theorems 3.3 and 3.9 to the rotor-router 
model, we are faced with two main problems. The first is to define an 
appropriate notion of convergence of integer-valued densities o"„ on 6n'^'^ to 
a real- valued density a on W^- The requirement that an take only integer 
values is of course imposed on us by the nature of the rotor-router model 
itself, since unlike the divisible sandpile, the rotor-router works with discrete, 
indivisible particles. The second problem is to find an appropriate analogue 
of Lemma ll. II for the rotor-router model. 

Although these two problems may seem unrelated, the single technique 
of smoothing neatly takes care of them both. To illustrate the basic idea, 
suppose we are given a domain A CM.'^, and let cr„ be the function on J^Z^^ 
taking the value 1 on odd lattice points in A'-, and the value 2 on even lattice 
points in A', while vanishing outside A'- We would like to find a sense in 
which cj„ converges to the real- valued density a = ^1a- One approach is 
to average an in a box whose side length L„ goes to zero more slowly than 
the lattice spacing: L„ J, while Ln/6n T oo as n t oo. The resulting 
"smoothed" version of an converges to a pointwise away from the boundary 
of ^. 

By smoothing the odometer function in the same way, we can obtain 



an approximate analogue of Lemma 1 . 1 for the rotor-router model. Rather 



than average in a box as described above, however, it is technically more 
convenient to average according to the distribution of a lazy random walk 
run for a fixed number a{n) of steps. Denote by (Xfc)fc>o the lazy random 
walk in 6n'^'^ which stays in place with probability ^ and moves to each of 
the 2d neighbors with probability Given a function / on iJ^Z*^, define its 
k-smoothing 

Skf{x)=E{f{Xk)\Xo = x). (53) 

From the Markov property we have SkSi = Sk+£- Also, the discrete Lapla- 
cian can be written as 

A = 26-^31 -So)- (54) 
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In particular, ASk = SfcA. 



4.1 Convergence of Odometers 

For n = 1,2,... let cj„ be an integer- valued function on 5„.Z'^ satisfying 
< (Tn < M. We assume as usual that there is a ball B CW^ containing the 
support of fT„ for all n. Let u be a function on supported in B satisfying 
( |27| ) and (28). In place of condition (30) we assume that there exist integers 
a{n) 1 oo with 5„a(n) | such that 

{Sa(n)<^n)°{x) ^ a{x), X ^ DC{a). (55) 

By the odometer function for rotor-router aggregation starting from ini- 
tial density o"„, we will mean the function 

Un{x) = 5,^ • number of particles emitted from x 

if (Tn{y) particles start at each site y. 

Theorem 4.1. Let Un he the odometer function for rotor-router aggregation 
on briJ^'^ starting from initial density an- If a, an satisfy (21)-(29), (31) and 
(55), then ^ s — 7 uniformly, where 

j{x) = -\x\'^ - g{x,y)a{y)dy 

and s is the least superharmonic majorant ofj. 

Given a function / on Sn'^'^, for a directed edge (x, y) write 

v/(.,,) = M_iw. 

On 

Given a function k on directed edges in 6n'^'^, write 

divK(x) = 7^^i^{x,y)- 



2d5n 



The discrete Laplacian on is then given by 

A fix) = div V/ = 5-2 ( ^ E f(y) - /(^) ) • 

The following "rescaled" version of \21\ Lemma 5.1] is proved in the same 
way. 
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Lemma 4.2. For a directed edge {x,y) in 6n'^'^, denote by K{x,y) the net 
number of crossings from x to y performed by particles during a sequence of 
rotor-router moves. Let 

u{x) = b\ ■ number of particles emitted from x during this sequence. 

Then 

Vu{x, y) = 6n{-2dK{x, y) + p(x, y)) (56) 
for some function p on directed edges in satisfying 

|p(x,y)| <4d-2 

for all directed edges {x,y). 

Proof. Writing N{x, y) for the number of particles routed from x to y, for 
any y, z ~ x we have 

\N{x,y)-N{x,z)\<l. 
Since u{x) = S^^^^^ N{x,y), we obtain 

5~'^u{x) -2d+l< 2dN{x, y) < 5^^m(x) +2d-l 

hence 

\Vu{x,y) + 2d6nti{x,y)\ = 5n\6n'^u{y) - 5-'^u{x) + 2dN{x,y) - 2dN{y,x)\ 

<{4.d-2)6n. □ 



Let Rn be the final set of occupied sites for rotor-router aggregation 
on J^Z'^ starting from initial density Since each site x starts with (T„(x) 
particles and ends with either one particle or none accordingly as x G Rn or 
X ^ Rm we have 

2d5nA\\ K = o-„ - 



Taking the divergence in ( 56 ) we obtain 

Au,„ = Ir^ -CTn + Sn div p. 

In particular, since < cr„ < M, and \p\ < 4d — 2, we have 

|Au„| <M + Ad 

on all of dr,'^'^. 



(57) 



(58) 
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Recall that the divisible sandpile odometer function has Laplacian 1 — cj„ 
inside the set Z)„ of fully occupied sites. The next lemma shows that the 
same is approximately true of the smoothed rotor-router odometer func- 
tion SkUn- Denote by and the probability and expectation operators 
for the lazy random walk {Xk)k>o defining the smoothing (53), started from 

Xq = X. 

Lemma 4.3. There is a constant Cq depending only on d, such that 

Co 



\ASkUn{x) - F^{Xk e Rn) + SkCTn{x)\ < 



Proof. Let k and p be defined as in Lemma 4.2, Using the fact that 
and A commute, we obtain from (57) 

ASkUnix) = P,(Xfc G Rn) - E,CT„(Xfc) + 6n div p(Xfc). (59) 



Since Vu and k are antisymmetric, p is antisymmetric by (56|. Thus the 
last term in (59) can be written 



Sn div p{Xk) = ^ ^ ^ P,(Xfc = y)p{y, z) 

y 2~y 



y zr^y 



= i E T^^^-^^k = y)- r.{X, = z))p{y, z) (60) 
y 2~?/ 

where the sums are taken over all pairs of neighboring sites y ^ z • 

We can couple lazy random walk started at x with lazy random walk 
X'f^ started at a neighbor x' of x so that the probability of not coupling in k 
steps is at most C/\/k -\- 1, where C is a constant depending only on d 
Since the total variation distance between the distributions of X^ and X'f, is 
at most the probability of not coupling, we obtain 



c{Xk 



z 



y zr^y 



)I=EE|IP- {Xk = y)-F^'{X', = y)\ 
y x'r^x 
^ 4dC 



Vk + l' 



Using the fact that \p{y,z)\ < Ad in (60), taking Cq = AdC completes the 
proof. □ 
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The next lemma shows that smoothing the rotor-router odometer func- 
tion does not introduce much extra error. 



Lemma 4.4. \SkUn 



Ur, 



<6iaMk + CQVk + l). 



Proof. From (54) we have 

\SkUn Ur, 



k-1 

- ^2 I'^i+i''^" 



Sj Un I 



j=0 

5: 



\ASjUn\. 



But by Lemma 4.3 



j=0 



\ASjUn\ < M + 



Co 



Summing over j yields the result. 
Let 

ln{x) = -\x\ 



□ 



GnSa{n)^n{x) 



and let s„ be the least superharmonic major ant of 7„. By Lemma the 
difference s„ — 7„ is the odometer function for the divisible sandpile on 
JnZ'^ starting from the smoothed initial density an = ^^(n)^^. Note that 
by Lemma 



3.4 



there is a ball O C M.'^ containing the supports of s — 7 and 
of Sn — 7n for all n. The next lemma compares the smoothed rotor-router 
odometer for the initial density (T.„ with the divisible sandpile odometer for 
the smoothed density fT„. 

Lemma 4.5. Let 17 C M*^ 6e a ball centered at the origin containing the 
supports of s — 'J and of Sn — 7n for all n. Then 

Sa{n)Un > Sn - In - Cqt'^ a{n)~'^ /"^ (61) 

on all of 5n'^'^, where r is the radius of Q.. 
Proof. Since 

A7n = -1 + Sa[n)(^n, 

by Lemma |4.3| the function 

f{x) = Sa(n)Un{x) + 7n(x) + Coa(n)"^/^(r^ - \xf) 

is superharmonic on Sn'^'^- Since / > 7n on il'', the function ipn = max(/, 7„) 
is superharmonic on Q''', hence -i/'n ^ s„ by Lemma 2.23 Thus / > on 
Q \ so (61) holds on rj'' and hence everywhere. □ 
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Lemma 4.6. Let be the set of occupied sites for rotor-router aggregation 
in SyiWj started froru initial density Cn- There is a ball ncR'^ with [JR° C 

n. 

Proof. By assumption there is a ball B C M.'^ containing the support of cr„ 
for all n. Let An be the set of occupied sites for rotor-router aggregation in 
(^nZ*^ started from initial density t{x) = MIx^b- From the abelian property 
we have Rn C An. By the inner bound of |2ll Theorem 1.1], if [26-'^jC{B)\ 
particles perform rotor-router aggregation starting at the origin in 5n 
the resulting set of occupied sites contains B'-; by the abelian property 
it follows that if we start with M [26n'^C{B)\ particles at the origin, the 
resulting set 0„ of fully occupied sites contains An. By the outer bound of 
|2H Theorem 1.1], 0,n is contained in a ball 0, of volume 2>MC{B). □ 

Let 

R^n^ = {x ^ Rn : y G Rn whenever ||x — < kSn}- 

The following lemma shows that the rotor-router odometer grows at most 
quadratically as we move away from the boundary of Rn. 

Lemma 4.7. For any k > 0, if x ^ Rn , then for any /3 > 

Un{x) < c{M + 4d)k'^6l 



where c = cp is the constant in Lemma 2.11 



Proof. If X ^ Rri \ then there is a point y Rn with |j; — y| < k6n. Since 



Un{y) = 0, we have by (|58[) and Lemma 2.17 



Unix) < c{M + 4d)\x - yl"^ . □ 



Proof of Theorem 4-1 By Lemma 3.4 there is a ball $7 containing the sup- 
port of s — 7 and of s„ — 7^ for all n. By Lemma 4.6 we can enlarge O if 
necessary to contain the support of Sa(_n)'^n for all n. By Lemma 4.3 the 
function 

(j){x) = Sa{n)Un{x) - Sn{x) + Jn{x) + Coa(n) "^^^ |xp . 

is subharmonic on the set 

Rn = {x £ Rn : y £ Rn whenever ||x — y||i < a{n)5n} 
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since {Xa(n) G Rn) = 1 for X G 

Given a point y ^ Rn, we have by Lemma 



4.7 



Uniy) < c{M + Ad)a{nf5l. 



Now by Lemma 4.4 



'5a(„)U„(y) < ( c(M + 4(i)a(n)2 + jMa(n) + Co^Ja{n) + 1 



The right side is at most Cia{n)'^5l, where Ci = c(M + Ad) + \M + Co^/l. 
Since Sn > 7n) we obtain for all y G J7 — i?„ 

(/-(y) < Cia(n)252 + Coa(n)-i/V2 (62) 



where r is the radius of 17. By the maximum principle, (62) holds for all 
y G 0,. Now from Lemma |4.5|we obtain 



(63) 



-Cor^a(n) < Sai^n)Un - Sn + 7n 



2.22 



and 



3.8 



we have 7^ — > 7 and — > s uniformly on Q. 

□ 



on all of (^nZ*^. 
By Lemmas 

Since a{n) j 00 and 5„a(n) | 0, we conclude from (63) that (»S'q(„)U„) 
s — 7 uniformly on 17. Since both Sa{n)'^n and s — 7 vanish outside $7, this 
convergence is uniform on W^. By Lemma 4.4 we have 



uniformly on ^riZ'^, and hence u,^ — > s — 7 uniformly on 
4.2 Convergence of Domains 



□ 



In addition to the assumptions of the previous section, in this section we 
require that 

For ah xeW^ either a{x) > 1 or a{x) = 0. 



We also assume that 



(64) 
(65) 



{a > 1} = {cr > ly. 
Moreover, we assume that for any e > there exists A^(e) such that 

If X G {cr > then a„(x) > 1 for all n > N{e); (66) 

and 

lix^{a> ly, then an{x) = for all n > N{e). (67) 
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Theorem 4.8. Let a and cr„ satisfy (^-@), @, ^ and (fg^-(fg7|j. 

For n > 1, let Rn he the final domain of occupied sites for rotor-router 
aggregation in 5n'^'^ started from initial density an- For any e > we have 
for all sufficiently large n 



(68) 



where ^ 

D = {s>-f}U{a> 1}°. 

Lemma 4.9. Fix e > and n > N{e/A). Given x i^*^/^ and 5n < p < e/4, 
let 

Mp{x) = #B{x,p)nRn. 

If Un < i^n"^ 0^ -S(x,p), t/ien 



M,{x)>(^l + ^^M,^sA^). 



Proof. By (65), we have {a > lY^^ C D^^^. The bah B(x,p) is disjoint 
from D*^/^, so by (67) we have cr„(y) = for all y G B{x,p). Since at least 
Mp-s,A^) particles must enter the ball B{x, p — we have 

un{y) > 6lNp-sSx). 

yedB{x,p-5n) 

There are at most Np{x) — J\fp-s„ix) nonzero terms in the sum on the left 
side, and each term is at most S^m, hence 

m{j\fp{x) - Mp-sAx)) > Mp^sA^)- n 



The following lemma can be seen as a weak analogue of Lemma 3.10 



for the divisible sandpile. In Lemma 4.11 below, we obtain a more exact 
analogue under slightly stronger hypotheses. 

Lemma 4.10. For any e > 0, if n is sufficiently large and x € Rn with 
X D''/'^ , then there is a point y G i?„ with \x — y\ < e/4 and 



Un{y) > 



e6n 



16 log {2u;d{e/46ny) 
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Proof. Take n > N{e/4) large enough so that 85„ < e. Let m be the 
maximum value of S~'^Un on ^(a;, e/4)''. Since x G Rn we have A/'o(x) = 1. 



Iteratively applying Lemma |4.9[ we obtain 

m 

For large enough n we have 



K/i{x) > f 1 + - 



J\fo{x) > exp 



16m6n 



M,/i{x) <#B{x,e/4y- <2ud 
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hence 



16m6n 

Solving for m yields the result. 



<log {2u;d{e/46n)''). 



□ 



Lemm a 4. 11. Fix e > and k > 4Cq, where Cq is the constant in 
Lemma 
isfies 



.3. There is a constant Ci such that if x G Rn, x D^^/^ sat- 

(69) 



and n is sufficiently large, then there exists y G 6n'^'^ with \x — y\ < ^ + 6n 
and 

SkUn{y) > SkUn{x) + —. 

Proof. Let Ci = c{M + Ad) + \M + Co\/2. Note first that by Lemma 
and (pi) 



un{x) > SkUn{x) - 5^ i^-^kM + CoVk + 1 

2r2 



4.4 



> c(M + 4d)r<5^ 



so a; G -R^f^ by Lemma 



4.7 



By Lemma |4.3[ we have 

ASkUn{y)>l-Skan{y)-Co/Vl^, y^P'n^. 

Note that Cq/V^ + 1 < 5- Take n > N{e/8) large enough so that M„ < e/8; 
then SkCTn vanishes on the ball B = i?(x,e/2)'' by (67). Thus the function 

f{z) = SkUn{z) - \\X - 
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~ (k) 

is subharmonic on i?„ := Kii DB, so it attains its maximum on the bound- 
ary. 



4.7 



If z S dR^^ , then by Lemma 
From Lemma 14.41 we obtain 

SkUn{z) < Cik'^51 < SkUn{x). 

Thus f{z) < f{x). Since x G Rn, it follows that / cannot attain its maxi- 

^ ^ tk) 

mum in RnUdRn at a point z G dRii , and hence must attain its maximum 
at a point y G dB. Since f{y) > f{x) we conclude that 

SkUniy) = f{y) + ^\x - 

>SkUn{x) + ^(^^y. □ 



Proof of Theorem 



Fix e > 0. Let u 



7 and D = {u > 0}. By 



Lemma 12.191 we have 

D, C D„ U {cj > 1}^^ 

for some r/ > 0. Since the closure of D^j is compact and contained in D, we 
have u > on for some m,j > 0. By Theorem 4.1 for sufficiently large 
n we have Un > u — ^m^j > on D^, hence C Rn- Likewise, by (66) we 
have {a > 1}^ C Rn for large enough n. Thus D]; C i2n- 

For the other inclusion, fix x G (5„Z^with x ^ D''. Since u vanishes 

for all 
< e/4. 
it follows that 



4.1 



we have n„ < ^ on 

r so that k6. 



on the ball B = B(x,e), by Theorem 
sufficiently large n. Let k = 4Co, and take n large enoug 

2 

Then SkUn < ^ on the ball 3e/4)''. By Lemma 
SkUn < Cik'^6n on the smaller ball B' = B{x,e/Ay'- 
y € B' we have 

SkUniy) ^ c.k'^/^+Hl 



4.11 



In particular, for 



Un{y) < 



^y{Xk = y) 

where C2 is a constant depending only on d. For sufficiently large n the right 
side is at most e(5„/161og (2a;(i(e/45„)'^), so we conclude from Lemma 
that X Rn- 



4.10 



□ 
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5 Internal DLA 



Our hypotheses for convergence of internal DLA domains are the same as 
those for the rotor-router model: o" is a bounded, nonnegative, compactly 
supported function on M"' that is continuous almost everywhere and > 1 
on its support; and {an}n>i is a sequence of uniformly bounded functions 
on ^nZ*^ with uniformly bounded supports, whose "smoothings" converge 
pointwise to a at all continuity points of a. 



Theorem 5.1. Let a and an satisfy (21)-(29), (31), (55) and (64)-(61). 



For n > 1 let In be the random domain of occupied sites for internal DLA 
in 6n'^'^ started from initial density an- If ^n log n I 0, then for all e > we 
have with probability one 



where 



D'^ G In CI D'^ ' for all sufficiently large n, 
D = {s>j}U{a> 1}°. 



(70) 



5.1 Inner Estimate 

Fix n > 1, and label the particles in J^Z*^ by the integers 1, . . . , m„, where 
nin = Z^xe^nZ'* '^nix). Let Xi be the starting location of the particle labeled 
i, so that 

#{i\xi = x} = an{x). 

For each i = 1, . . . ,m„ let {Xl)t>o be a simple random walk in (5„Z'^ such 
that Xq = Xi, with X* and independent for i ^ j. 
For z G (JnZ'^ and e > 0, consider the stopping times 

Ti=mf{t>0\Xl = z}; 
Ti=mi{t>0\Xl^D-^}; 

u^ = mi{t>0\Xl^{Xi^yr},}. (71) 

The stopping time v"^ is defined inductively in i with = 0. We think 
of building up the internal DLA cluster one site at a time, by letting the 
particle labeled i walk until it exits the set of sites already occupied by 
particles with smaller labels. Thus z^* is the number of steps taken by the 
particle labeled i, and X^^ is the location where it stops. 
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Fix z G U'l and consider the random variables 




i=l 



These sums can be interpreted in terms of the following two-stage procedure. 
During the first stage, we let each particle walk until it either reaches an 
unoccupied site or exits D];. In the second stage, we let each particle that 
did not already exit D'^ continue walking until it exits D];. Then counts 
the number of particles that visit z during both stages, while counts the 
number of particles that visit z during the second stage. In particular, if 

< Ai^, then z was visited during the first stage and hence belongs to the 
occupied cluster /„. 

The sum is difficult to estimate directly because the indicator random 
variables in the sum are not independent. Following [20], we can bound 
by a sum of independent indicators as follows. For each site y S D];, let 
{Y^^)t>o be a simple random walk in 6n'^'^ such that Yq = y, with and Y^ 
independent for x / y. Let 



Thus counts the number of walks Y^ that hit z before exiting D'^ . Since 
the sites X^^ are distinct, we can couple the walks {Y^} and {X'^} so that 




where 



■y 



inf {t >0\Yy = z}; 
inf {t >0\Yy ^D'-]. 



Define 



fn,e{z) = gn,e{z, z) E {M^ - L,), 



where 



5n,.(y, z) = E#{t < ry I Y^ = z] 



is the Green's function for simple random walk in 6n'^'^ 



stopped on exiting 
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D'-. Then 



\i=l y&D]} 

= gn,.{z,z) ^ K(y) - 1) P (r| < ri') 

= ^K(y)-i)5n,.(y,^) (72) 

where in the last step we have used the identity 

P (rv < tV) = ^"'^i^' ■ (73) 

9n,e{z,z) 

Thus fn,e solves the Dirichlet problem 

A/„,, = 5-2(1 -a„), onD^- (74) 

/n,. = 0, ona^^ (75) 

Note that the divisible sandpile odometer function Un for the initial 
density (T„ solves exactly the same Dirichlet problem, with the domain D'^ 



in (74) and (75) replaced by the domain D„ = {un > 0} of fully occupied 
sites. Our strategy will be first to argue that since the domains Dn and D are 
close, the solutions to the Dirichlet problems in Dn and D]; should be close; 



next, by Theorem 3.3 since Un ^ u := s — j it follows that the functions 
fn,e and u are close. Since u is strictly positive in D^, we obtain in this way 
a lower bound on /„^e, and hence a lower bound on E(A4e — i^). Finally, 
using large deviations for sums of independent indicators, we conclude that 
with high probability < Me, and hence with high probability every point 
z £ D'-^ belongs to the occupied cluster I„. The core of the argument is 



Lemma 5.4 below, which gives the desired lower bound on fn,e 



In order to apply Theorem 3.3 we must have discrete densities which 
converge pointwise to a at all continuity points of a. Recall, however, that in 
order to allow a to assume non-integer values, we have chosen not to assume 
that (T„ converges to cr, but rather only that the smoothed density Sf^(^n)'^n 



converges to a; see (55). The next lemma shows that this initial smoothing 



step does not change the divisible sandpile odometer by very much. 

Lemma 5.2. Let k >0 be an integer, and let Un (resp. Un) be the odometer 
function for the divisible sandpile on 5n'^'^ started from initial density an 
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(resp. SkCTn)- If CTn has finite support and < cj„ < M, then 



1 



un\ < -kM6f^. 
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Proof. From (54) we have the identity Sif = / + -f^f- By induction on A;, 



it follows that 
where 



Wk 



j=0 



Since < 5j(T„ < M for all j, we have < Wfc < ^kMS"^. 

Let t'n (resp. P„) be the final mass density on (J^Z^^ for the divisible 
sandpile started from initial density an (resp. SkCTn)- Then 

i^n = SkCTn + Au„ 

= cj„ + A(n„ + Wk) . 



Since t'n < 1 and m„, + Wfc ^ 0, Lemma 3.2 yields Un < 'Un + u^fc- Likewise 
= S'feO-n + A (un -Wk + ^kM6l 



Since m„ < 1 and Un — Wk + \kM5n > 0, Lemma 3.2 yields Un < Un — Wk + 
\kM6l. □ 



The next lemma adapts Theorem 3.3 to our current setting. The hy- 
potheses are identical to those of Theorem 3.3 except that the convergence 
of smoothed densities (55) replaces (30). 

Lemma 5.3. Let Un be the odometer function for the divisible sandpile on 
with initial density an, and let Dn = {un > 0}. If a and an satisfy 

s — 7 uniformly. Moreover, for 



(2^-ll2S^, and ([£5p, then ^ u := 
any e > we have Dn C -D'' for all sufficiently large n. 

Proof. Let n„ be the odometer function for the divisible sandpile on 
started from initial density Sa[n)'^n- Then ^ u uniformly by Theo- 
rem 3.3 By Lemma 5.2 we have for all z G M'^ 



n°(z)-ti^(z) <^-Ma[n)5l 



1 
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Since the right side tends to zero as n | oo, we obtain u uniformly. 

Let /3 > be the minimum value of n on Z^^. Taking n large enough so 
that \un — u''\ < (3/2, we have n„ > (3/2 on D];, hence Dn C -D''. □ 

Lemma 5.4. Fix e > 0, and let P > be the minimum value of u = s — j 
on D^. There exists < rj < e such that for all sufficiently large n 

fn,r,{z)>\(35-\ ZGDl 

Proof. Since u is uniformly continuous on D and vanishes outside D, we 
can choose rj > small enough so that u < (3/6 outside Drj. Since u > (3 
on D^, we have r] < e. Let u„ be the odometer function for the divisible 
sandpile on started from initial density c^i, and let Dn — ^ 0} 

be the resulting domain of fully occupied sites in br^'^- We have u 
uniformly by Lemma 5.3 so for n sufficiently large we have \un — u'''\ < /3/6, 
hence n„ < /3/3 on 91)^. 



By Lemma 5.3 we have C Dn for all sufficiently large n. Now (74) 
implies that the difference S'^fn,ri — Un is harmonic on Z)^, so it attains its 
minimum on the boundary. Since fn,r] vanishes on 51?^, we obtain for z E D':/^ 

blfn,r,{z) > Un{z) - - 

> uiz) . 

Since u > (3 in D'^, we conclude that 5nfn,r) > (3/2 in D'^. □ 
Lemma 5.5. We have 

ELr - ^ 



9n,ri^Z, z) 

Moreover 

EA^^ < M- ^ 



9n,ri{z, z) 



Proof. By (73) and the symmetry of gn,rp we have 

EL, = P(r| < ry) (76) 

^ ^ 9n,r,{V,z) 
y&D- 9n,v{z,z) 



9n,rj{z,Z, ^^^^ 
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The sum in the last hne is KzTqd-^. 

To prove the inequahty for ETW^, note that EA^^ is bounded above by 
M times the sum on the right side of (76). □ 

The next lemma, using a martingale argument to compute the expected 
time for simple random walk to exit a ball, is well known. 

Lemma 5.6. Fix r > and let B = B{o,ry'' C 6n'^'^, and let T be the first 
hitting time ofdB. Then 

Proof. Since 5^^|Xtp — t is a martingale with bounded increments, and 
EoT < cxD, by optional stopping we have 

EoT = 6-^Eo\Xt\^ = 6-\r + 0{5n)f. □ 



The next lemma, which bounds the expected number of times simple 
random walk returns to the origin before reaching distance r, is also well 
known. 

Lemma 5.7. Fix r > and let B = B{o,rY'' C 5jtTj'^, and let Gb be 
the Green's function for simple random walk stopped on exiting B. If n is 
sufficiently large, then 

r 

Gb{o,o) < log— . 

On 

Proof. In dimension d >3 the result is trivial since simple random walk on 
6n'^'^ is transient. In dimension two, consider the function 

f{x) = Gb{o,x) - gn{o,x) 



where gn is the rescaled potential kernel defined in (26). Since / is harmonic 
in B it attains its maximum on the boundary, hence by Lemma 2.21 we have 
for X £ B 

2 f6l' 
f{x) < — logr + O 



Since gn{o, a) = ^ logf^m the result follows on setting x = a. □ 
We will use the following large deviation bound; see [21 Cor. A. 14]. 
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Lemma 5.8. If N is a sum of finitely many independent indicator random 
variables, then for all X > 

F{\N - EN\ > XEN) < 26-"^^^ 

where cx > is a constant depending only on A. 

Let 

in = {Xl^W' < f'} C In 



where is given by (71), and 

f' = inf > 0\Xi in-']. 

The inner estimate of Theorem |5.1| fohows immediately from the lemma 
below. Although for the inner estimate it suffices to prove Lemma 5.9 with 
/„ in place of /„, we will make use of the stronger statement with /„ in the 
proof of the outer estimate in the next section. 

Lemma 5.9. For any e > 0, 

¥{^D'^ C In for all but finitely many n) = 1. 

Proof. For z € D'^, let Szi^) be the event that z ^ In- By Borel-Cantelli it 
suffices to show that 

F{£z{n)) < oo. (77) 



By Lemma 



2.19 



since D = D U {a > 1}° we have 
A C A' U {o- > 1},/ 



for some e' > 0. By (66), for n > N(e') the terms in (77) with z G {cr > 1}^ 
vanish, so it suffices to show 



EE 

n>l zS-D-', 



< oo. 



(78) 



By Lemma 5.4 there exists < rj < e' such that 



1 



z £ D -j, 



(79) 
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for all sufficiently large n, where /3 > is the minimum value of u on 
Fixing z £ since < we have 



:(n)) < F{M:r, = Lr,) 

< F{M^r, < a) + ¥{Lr, > a) 



(80) 



5.8 



for a real number a to be chosen below. By Lemma 
are sums of independent indicators, we have 

F {Lrj > (1 + A) E Lrj) < 2e~^^'^^'' 
P {M^ < (1 - A) E Mr^) < 26"^^^^" 



since and A^^ 



(81) 



where cx depends only on A, chosen below. Now since z E D^' and D is 
bounded by Lemma 2.15 we have 

B{z,e' -rj) C Dr, C B{z,R) 



hence by Lemma 5.6 



1 fe'-r] 
2 



< Er^^ < 



(82) 



By Lemma 5.5 it follows that KAirj < MR'^ /6'^gn,ri{z, z). Taking 

(3 



a = EL„ + 



'^^?i9n,rj{z,z) 



in (|80|), and letting A = f3/4MR^ in {jSl), we have 



(83) 



hence 



a > (1 + A) EL^. 



Moreover, from (79) we have 



hence by ( 83 1 



a < EM 



^5l9n,r,{z,z] 



< (1 - X)EM^. 
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Thus we obtain from (80) and (|8lj) 



By Lemmas 5.5 and 5.7 along with (82) 



EL. 



> 



1 fe' -rj 



2\ 5n J log(i?/<5. 



Using (84), the sum in (78) is thus bounded by 



n>l zGD-'; 



n>l 



cx{e' - 7?)^ 
'25llog{R/6r, 



(84) 



< oo. □ 



5.2 Outer Estimate 

For X G Z'^ write 

Q{x,h)={yeZ'^ : \\x - y\\^ < h} 

for the cube of side length 2h + 1 centered at x. According to the next 
lemma, if we start a simple random walk at distance h from a hyperplane 
H C Ij'^, then the walk is fairly "spread out" by the time it hits H, in the 
sense that its chance of first hitting H at any particular point z has the 
same order of magnitude as z ranges over a (d — 1) -dimensional cube of side 
length order h. 

Lemma 5.10. Fix y ^TJ^ with yi = h, and let T be the first hitting time of 
the hyperplane H = {x £ 'L'^\xi = 0}. Let F = H n Q{y, h). For any z £ F 
we have 

FyiXr = z)> ah^-'^ 
for a constant a depending only on d. 
Proof. For fixed w £ H, the function 

fix) = F,{Xt = w) 

is harmonic in the ball B = B{y, h). By the Harnack inequality |19| Theorem 
1.7.2] we have 

f{x)>cf{y), xeB{y,h/2) 
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for a constant c depending only on d. By translation invariance, it follows 
that for w' £ H f] B{w, h/2), letting x = y + w — w' we have 

FyiXr = w') = fix) > cfiy) = cFy {Xt = w). 

Iterating, we obtain for any w' £ H f] Q{w, h) 



{Xt = w')>c 



{Xt = w) 



(85) 



Let T' be the first exit time of the cube Q{y,h — 1). Since F is a boundary 
face of this cube, we have {Xt' G F} C {Xt € F}. Let zq be the closest 
point to y in H. Taking w' = zq, we have w' S Hr]Q{w, h) whenever w € F. 
Summing (85) over w £ F, we obtain 

{2h + if-^c-^Fy {Xt = zo) > Fy{XT G F) > F{Xt' ^ F) = ^- 



Now for any z £ F, taking w = zq and w' = z in (85), we conclude that 



{Xt = z)>c^^ Fy {Xt = zo) > 



Id 



-(2/1 + 1 



□ 



Lemma 5.11. Let h,p he positive integers. Let N he the numher of parti- 
cles that ever visit the cuhe Q{o, p) during the internal DLA process, if one 
particle starts at each site y — Q{o,p + h), and k additional particles 
start at sites yi, . . . ,yk ^ Q(o, p + h). If k < \h'^, then 

N < Binom(/c,p), 

where p < 1 is a constant depending only on d. 

Proof. Let Fj be a half-space defining a face of the cube Q = Q{o,p), such 
that dist{yj,Fj) > h. Let Zj G Fj he the closest point to yj in Fj, and let 
Bj = Q{zj, h/2) n Fj. Let Zj be the random set of sites where the particles 
starting at yi, ... ,yj stop. Since #Bj > ^h'^ > 2k, there is a hyperplane Hj 
parallel to Fj and intersecting Bj, such that 

#Bj n Hj n Zj^i <\h'^^\ (86) 

Denote by Aj the event that the particle starting at yj ever visits Q. On this 
event, the particle must pass through Hj at a site which was already occupied 
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by an earlier particle. Thus if {Xt}t>o is the random walk performed by the 
particle starting at yj, then 



where T is the first hitting time of Hj. By Lemma 5.10 every site z G BjHHj 
satisfies 

F{Xt = z)> ah^-'^. 

From ([86|, since #5^ n Hj > /i"'"^ we obtain F{Xt Zj^i) > a/2, hence 

P(Aj-|.^,_i) <p 

where p = 1 — a/2, and .T^j is the a- algebra generated by the walks performed 
by the particles starting at yi, . . . , y^. Thus we can couple the indicators 
with i.i.d. indicators Ij > of mean p to obtain 

k k 
TV = ^ 1a^, < ^ Ij = Binom(A;,p). □ 



The next lemma shows that if few enough particles start outside a cube 
Q{o,3p), it is highly unlikely that any of them will reach the smaller cube 
Q{o,p). 

Lemma 5.12. Let p,k be positive integers with 

A:<i(l-//2.)V 



where p < 1 is the constant in Lemma 5.11. Let N be the number of parti- 
cles that ever visit the cube Q{o, p) during the internal DLA process, if one 
particle starts at each site y G X'^ — Q{o,3p), and k additional particles start 
at sites yi, ■ ■ ■ ,yk ^ Q{o, 3p). Then 

F{N > 0) < coe-'^'P 

where co,ci > are constants depending only on d. 

Proof. Let Nj be the number of particles that ever visit the cube Qj = 
Q{o, Pj), where 

p, = {2 + p^/''^)p. 
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Let kj = p'^'^k, and let Aj be the event that Nj < kj. Taking h = pj — pj+i 
in Lemma |5.11 



since 



we obtain 
Hence 



k,<\j^/'{l-p'/'Yp'=\h' 



iVj+iU, < Bmom{Nj,p). 



F{Aj+i\Aj) > F{Bmom{k j,p) < kj+i) 
> 1 - 2e-^^^ 



(87) 



where in the second hne we have used Lemma 5.8 with A = ^/p — p. 
Now let 

. ^ ^ logfc - log/) 
L log(l/p) 

so that p^/'^p < kj < p. On the event Aj, at most p particles visit the 
cube Q{o,2p). Since the first particle to visit each cube Q{o,2p — i) stops 
there, at most p — i particles visit Q(o,2p — i). Taking i = p we obtain 
F{N = 0) > F{Aj). From ([87]) we conclude that 

P(iV = 0) > P(^i) P (yl2|^i) • • • IP (AjlAj^i) > 1 - 2je-''P. 

The right side is at least 1 — coe"'^^'' for suitable constants cq, ci. □ 



Proof of Theorem 5.1. The inner estimate is immediate from Lemma 5.9 
For the outer estimate, let 



#{l < i < mn\y' > f'} 



5.11 



For e > 



be the number of particles that leave D ' before aggregatin g to t he cluster. 

Let Kq = ^ (^ '^^^3^ — ~^ ' where p is the constant in Lemma 
and no > 1, consider the event 

Fno = {Nn < Ko{e/6ny for all n > no}. 



sition 2.12 



By (64) and (65), the closure of D contains the support of a, so by Propo- 

d J2 <Tn{x)^ [ a{x)dx = C{D). 
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Moreover, by Proposition 2.12[ i), for sufficiently small rj we have 
Taking n large enough so that 6n < rj, we obtain 



Thus for sufficiently large n, on the event D]^ C In we have 

Nn<mn- 



1 



<mn- 6-''C[D) + -Koe'^d, 



where in the last line we have used (88). From Lemma |5.9| we obtain 
P {Fno) > 1P(5-^ ^ In for ah n > no) T 1 



(89) 



as no T oo. 

By compactness, we can find finitely many cubes Qi, . . . ,Qj of side 
length £ = 2e/3\/d centered at points xi,...,Xj S d{D'^), with d{D'^) C 
UQi. Taking k = [Ko{e/6nf\ and p = \e/6 n\ in Lemma 5.121 since the 
cube of side length 3ii centered at Xi is disjoint from D, we obtain for n > uq 

r{{Q^ n /„ / 0} n Fn,) < 

Since {in ^ D'^''-^ C Ui=i {<3i ^ -^n 7^ O}, summing over i yields 

By Borel-Cantelli, if G is the event that ^ D'^'' for infinitely many n, 
then P(F„,(, n G) = 0. From ^ we conclude that P(G) = 0. □ 



6 Multiple Point Sources 



This section is devoted to proving Theorem 1.4 We further show in Propo- 



sition that the smash sum of balls arising in Theorem 1.4 is a classical 
quadrature domain. In particular, the boundary of the smash sum of k disks 
in lies on an algebraic curve of degree at most 2k. 
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6.1 Associativity and Continuity of the Smash Sum 

In this section we establish a few basic properties of the smash sum Q. 
In Lemma 16.11 we show that the smash sum is associative. In Lemma 16.21 
we show that the smash sum is continuous with respect to the metric on 
bounded open sets A, B CW^ given by 

dA{A,B) = £{A A B) 

where A A B denotes the symmetric difference A D B — A f] B, and C is 



Lebesgue measure on M*^. Lastly, in Lemma 6.4 we show that the smash 



sum is also continuous in a variant of the Hausdorff metric. 

Lemma 6.1. Let A,B,C C M*^ be bounded open sets whose boundaries have 
Lebesgue measure zero. Then 

{A®B)®C = A®{B®C) 

= AU BUCU D 



where D is given by (21) with o" = 1^ + 1^ + Ic. 
Proof. Let 

j{x) = -\x\^ -G{1a + 1b){x) 

and 

7(2;) = -|x|2-G(lAeiJ + lc)(a:). 

Let u = s— 7 and u = s— 7, where s, s are the least superharmonic majorants 
of 7, 7. Then 

{A® B) ® C = {A® B) U C U {u > 0} 

= A U 5 U {m > 0} U C U {u > 0} 

= AuB\JCU{u + u> 0}. (90) 

Let Un be the final mass density for the divisible sandpile in (5„Z'^ started 
from initial density 1^- + Is-, and let be the corresponding odometer 



function. By Theorem 3.3 we have Un ^ u as n ^ 00. Moreover, by 
Theorem 13.91 we have 

i^°{x) ^ lA^Bix) (91) 
for all X ^ d{A®B). Let m„ be the odometer function for the divisible sand- 



pile on Sji"^ started from initial density Vn + Ic--- By Proposition 2.12 'i) 



the right side of (91) is continuous almost everywhere, so by Theorem 3.3 
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we have ^ -u as n — > oo. On the other hand, by the abehan property, 
Lemma 13. 11 the sum Un + Un is the odometer function for the divisible sand- 



pile on 6n'^'^ started from initial density 1^- + Ib-- + lc--i so by Theorem 
we have Un + Un ^ u '■= s — j, where 

j{x) = -|x|2-G(1a + 1b + 1c)(x) 

and s is the least superharmonic majorant of 7. In particular, 

u + u = lim (Un + Un) = u, 



3.3 



so the right side of (90 ) is equal to AU B U C U D. 



□ 



Lemma 6.2. Let Ai, A2, B be bounded open subsets ofW^ whose boundaries 
have measure zero. Then 

dA{Ai®B,A2®B) <dA(Ai,^2)- 



Proof. Let A = AiU A2, and A[ = A- Ai for i = I, 2. Then by Lemma [(U 
Ae B = {A'i® Ai) e B = A'i® {Ai © B). 



By Corollary 2.13 



C{A ®B) = C{A'i) + C{A^ © B) 



hence 



dA{AeB,AieB) = C{A'i). 
By the triangle inequality, we conclude that 

dA{Ai(BB,A2®B) <C{A[)+C{A'2) = dA{Ai,A2). □ 



The following lemma shows that the smash sum of two sets with a small 
intersection cannot extend very far beyond their union. As usual, A^ denotes 
the outer e-neighborhood ([8| of a set A C M'^. 

Lemma 6.3. Let A, B C be bounded open sets whose boundaries have 
measure zero, and let p = {C{Ari B))^/^ . There is a constant c, independent 
of A and B, such that 

A®BC {AUB)''P. 
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Proof. Let 7,5 be given by ([I9j) and (20) with a = 1^1 + 1^, and write 
u = s - J. Fix X e {A ® B) - {AU B) and let r = dist(x,yl U B). Let 



B = B{x,r/2). By Lemma |2.2[ iii), s is harmonic in A (B B — {AU B) , so hy 
Lemma 12.41 the function 

w{y) = u{y) - \x- ?/|^ 

= s[y) + + G{Ia + lB)[y) -\x- y\^ 

is harmonic on the intersection {A® B) r\ B; hence it attains its maximum 
on the boundary. Since w{x) > the maximum cannot be attained on 
d{A (B B), so it is attained at some point y £ dB, and 



u{y) > w{y) + ^ > ^• 



(92) 



If z is any point outside A® B, then by Lemma 2.6 and Lemma |2.18 
with A = 4(i, there is a constant c' > 1 such that u < c'h? on B{z, h) for all 



h. Taking h = r/2vc', we conclude from (92) that B{y, h) <Z A® B. Since 



B{y, h) is disjoint from j4 U i?, we have by Corollary 2.13 

C{A UB) + LOdh'^ < C{A ®B) = C{A) + C{B), 
hence tOdh'^ < C{An B). Taking c = 2-v/c' I'^J'^ yields r = cio^J'^h < cp. □ 



The next lemma, together with Lemma 2.8 shows that the smash sum 
is continuous with respect to the metric 

dio{A, B) = ini {r > \ Ar C B C A" and Br C A C B''} 

on bounded open sets A,Bc M"*. Here A^ and A"' denote the inner and 
outer neighborhoods ([8| of A. This metric, which we call the inner-outer 
Hausdorjf metric, is related to the usual Hausdorff metric dn by 

dioiA, B) = max {dn {A, B) , dniA'^, B')) . 

Lemma 6.4. Let A,B C be bounded open sets whose boundaries have 
measure zero. For any e > there exists rj > such that 



{A © B), C Ar,® Br^ C A"^ ® Bi C {A® Bf 



(93) 



Proof. By Lemma 2.19 since A® B = AVJ B \J D, we have 

{A ® B)e C A,> U B^> U A' 
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for some e' > 0. By Lemma 2.16 m) with a = 1a+^b and (t„ = 1aj^/„+Ibi/„, 
for sufficiently small rj we have C ylr,©i?^. This proves the first inclusion 
in (|93j). 

The second inclusion is immediate from Lemma 12.81 
For the final inclusion, write A' = — A and B' = — B. Since 
C{A') i C{dA) = and £{B') i C{dB) = as t/ j 0, for small enough r] we 



have by Lemmas 6.1 and 6.3 

A"^ (BB"^ = AeBeA' eB' 

C {{AeB)uA'^UB'^Y 
C{A®BY. □ 



6.2 Smash Sums of Balls 



In this section we deduce Theorem 1.4 from our other results, using the 



main results of ^20j and |21j. The following theorem collects the results 
of |2H Theorem 1.3] for the divisible sandpile, |2H Theorem 1.1] for the 
rotor-router model, and |20| Theorem 1] for internal DLA. 

Theorem 6.5. Fix A > and a sequence 6n I 0, and let Dn be the domain 
of fully occupied sites for the divisible sandpile in b.^^L'^ started from mass 
m = [\5-'^\ at the origin. Likewise, let Rn,In be the domains of occupied 
sites in b^^'^ formed respectively from the rotor-router model (with initial 
rotors configured arbitrarily) and internal DLA, starting with m particles at 
the origin. For any e > 0, we have 

B'^ C Dn, Rn C B'' ' for all sufficiently large n, 

where B is the ball of volume A centered at the origin in M'^. Moreover, if 
bn < 1/n, then with probability one 

B'^ C In C B^" for all sufficiently large n. 



Proof of Theorem I.4 For i = 1, . . . , A; let DJ^, R^, 1^ be the domain of oc- 
cupied sites in (5„Z'^ starting from a single point source of |_(5~'^AjJ particles 
at x'^. By Theorem 16.51 for any r/ > we have with probability one 



{BiY^ C Dn, Rn C (Sj)''-- for all i and ah sufficiently large n; (94) 
and 'd bn < l/?^, then with probability one 

(-Bj)5j C InC {Biy for all i and all sufficiently large n. (95) 
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Next we argue that the domains Dn, Rn, In can be understood as smash 
sums of Dl^,K''^,I^ as i ranges over the integers l,...,k. By the abehan 
property |7], the domain /„ is the Diaconis-Fulton smash sum in dnZ"^ of 
the domains Likewise, if ri is an arbitrary rotor configuration on Sn'^'^, 
let be the smash sum of i?^ and formed using rotor-router dynamics 
with initial rotor configuration ri, and let r2 be the resulting final rotor 
configuration. For i > 3 define 5"^ inductively as the smash sum of SJ^^ 
and -R^ formed using rotor-router dynamics with initial rotor configuration 
rj_i, and let ri be the resulting final rotor configuration. Then Rn = S^. 



Finally, by Lemma 3.1 the domain Dn contains the smash sum of domains 
formed using divisible sandpile dynamics, and Dn is contained in the 
smash sum of the domains Dt U DDI, . 



Fixing e > 0, by Theorem 1.3 and Lemma 2.8, it follows from (94) and 



( 95 ) that for all sufficiently large n 



c 



(96) 



where A is the smash sum of the balls {Bi)n, and A is the smash sum of the 



balls (-Bj)''. By Lemma 6.4 we can take rj sufficiently small so that 

A/2 CACAC D''^ 

where D 



Bk ■ Together with ( 96 ) , this completes the proof. □ 



Proposition 6.6. Fix xi, . . . ,Xk E M*^ and Ai, . . . , > 0. Let Bi he the 

ball of volume Aj centered at Xi . Then 



/ h{x)dx < 'S~^ \ih{xi) 



for all integrable superharmonic functions h on Bi 



Bk- 



We could deduce Proposition |6.6| from Proposition 2.11 by integrating 
smooth approximations to h against smooth densities an converging to the 
sum of the indicators of the B^. However, it more convenient to use the fol- 
lowing result of M. Sakai, which is proved by a similar type of approximation 
argument. 

Theorem 6.7. ^ Thm. 7.5] Let CW^ he a hounded open set, and let 
a he a hounded function on supported on Cl satisfying a > 1 on Q.. Let 



'~i,s,D he given hy (19)-(21). 



h{x)dx < / h{x)a{x)dx 
D Jd 
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for all integrable superharmonic functions h on D. 

We remark that the form of the obstacle in p7j is superficially different 
from ours: taking w = 1 and uj = a \x\ section 7 of [27^, the obstacle is given 
by 

V'(x) =GlB-Ga 



for a large ball B, rather than our choice (19) of 

'2 Ga. 



7(x) 



Note, however that 7 — -0 is constant on B by (17) and (18). Thus if B is 



sufficiently large, the two obstacle problems have the same noncoincidence 



set D by Lemmas 2.9 and 2.15 



To prove Proposition 6.6 using Theorem 6.7 we must produce an appro- 
priate density a on strictly exceeding 1 on its support. We will take a 
to be twice the sum of indicators of balls of half the volume of the Bj . 



Proof of Proposition 6.6 Let B'- be the ball of volume Aj/2 centered at Xi 
and consider the sum of indicators 



a = 2^1^,. 



i=l 



Let j^SjD be given by (19)-(21). By Theorem 6.7 and the mean value 



proprerty for superharmonic functions (12), we have 

k 



/ h{x)dx < 2 ^ / h{x)dx < ^ Ai/i(; 

1 = 1 -^^'r i=l 



for all integrable superharmonic functions h on D. 

It remains to show that D = i?i © . . . © i^^. By Lemma 2.7 we have 
B[ d D for all i, hence by Lemma 6.1 



D = B[® B[ 



B'k®Bl 



By Lemma 2.14 we have Bi = B'^ (B B'^, completing the proof. 



□ 



From Proposition 6.6 it follows |13| sec. 6] (see also \14:\ Lemma 1.1(a)]) 
that if xi, . . . , Xfc are distinct points in M^, and Bi is a disk centered at Xi, 
the boundary of the smash sum Si © ... © lies on an algebraic curve of 
degree 2k. The methods of [13] rely heavily on complex analysis, which is 
why they apply only in dimension two. Indeed, in higher dimensions it is 
not known whether classical quadrature domains have boundaries given by 
algebraic surfaces p9|, Ch. 2]. 
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7 Potential Theory Proofs 



In this section we collect the proofs of the background results stated in 
section [2] For more background on potential theory in general, we refer the 
reader to [31 18] ; for the obstacle problem in particular, see IH [10] . 

7.1 Least Superharmonic Majorant 

Proof of Lemma \2.S\ (i) Let / > 7 be continuous and superharmonic. Then 
f > s. By the mean value property (12), we have 

f>Arf> ArS. 

Taking the infimum over / on the left side, we conclude that s > ArS. 
It remains to show that s is lower-semicontinuous. Let 

uj{j,r)= sup - liy)]. 

^ ,\x—y\<r 



Since 

ArS > Arj > 7 — '^(7, r) 

the function Aj-s + 1^(7, r) is continuous, superharmonic, and lies above 7, 
so 

ArS < s < ArS + a;(7, r). 

Since 7 is uniformly continuous, we have u;(7, r) | as r | 0, hence ArS — > s 
as r I 0. Moreover if rg < ri, then by Lemma |2.1[iii) 



Ar^S = lim ArgArS > lim Ar-^ArS = Arj^S. 
r^O r^O 

Thus s is an increasing limit of continuous functions and hence lower- 
semicontinuous . 

(ii) Since s is defined as an infimum of continuous functions, it is also 
upper-semicontinuous . 

(iii) Given x £ D, write e = s{x) — 7(x). Choose 6 small enough so that 
for ally e B = B{x, 5) 

|7(x) -7(y)| < I and \s{x) - s{y)\ < ^. 

Let / be the continuous function which is harmonic in B and agrees with 
s outside B. By Lemma |2.1[ v), / is superharmonic. By Lemma |2.1[ i), / 
attains its minimum in i? at a point z G dB, hence for y B 

f{y) > f{z) = s{z) > s{x) - ^ = 7(x) + ^ > 7(y). 
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It follows that / > 7 everywhere, hence f > s. From Lemma |2.1[ ii) we 
conclude that f = s, and hence s is harmonic at x. □ 

Proof of Lemma \2.4\ Suppose B{x,r) C il. Since for any fixed y the func- 
tion f{x) = g{x,y) is superharmonic in x, we have 

Ga{x) = / a{y)g{x,y)dy 

> [ a{y)Arf{y)dy 
= ArGa{x). □ 



Proof of Lemma \2.5\ (i) By Lemma 2.4 since a is nonnegative, the function 
7(x) + |xp = —Ga{x) is subharmonic on M'^. 

(ii) Let B = B{o, R) be a ball containing the support of a. By (17) and 
(|18|), for X B we have 



= CdR^ - Glsix) 

where C2 = 1 — 2 log R and Cd = d > 3. Hence for x & B we have 

7(x) - (M - l)\x\'^ = -Ga{x) - M|x|2 

= G{M1b - a){x) - CdMR^. 



Since a < MIb on by Lemma 2.4 the function 7 — (M — l)|xp is super- 
harmonic in B DQ. Since this holds for all sufficiently large R, it follows 
that 7 — (Af — l)|xp is superharmonic on all of □ 



Proof of Lemma \2.()\ (i) By Lemmas |2.2[ i) and |2.5[ i), the function 

u — = s — (7 -|- \x\'^) 

is the difference of a superharmonic and a subharmonic function, hence 
superharmonic on M'^. 

(ii) With Ar defined by (12), we have 



1 



JB{o,r) 



|xp + 2x ■ y + |y|^) dy 



\x\ + 



1 



{diVdt'^-^)t^ dt 



\x? + 



^dr" Jo 
dr'^ 



d + 2 
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By Lemma 2.5 we have 



ArS + —— > Ar-f + Ar|xp 



\x\ > 7. 



Since ArS is continuous and superharmonic, it follows that ArS + > s, 
hence 



d + 2 



> s + \x\ 



Thus s + |xp is subharmonic on W^, and hence by Lemma 2.5 ii) the function 
u + M|xp = (s + - (7 - (M - 



is subharmonic on J7. 



□ 



Proof of Lemma 2.7 If cr > 1 in a ball B = B{x,r), by (17) and (18), for 
y & B we have 



7(y) 



Ga(y) = -Cdr2-G(a-lB)(y). 



By Lemma |2.4| it follows that 7 is subharmonic in B. In particular, s > 7 
in B, so X £ D. □ 



Proof of Lemma 2.8. Let 



S = S2 + G(cr2 - (Tl). 



Then s is continuous and superharmonic, and since Si{x) > — |x| 
we have 

six) > — |xp — Gai{x) 

hence s > si. Now 

U2 — Ul = S2 — Sl + G((T2 — CJl) = S — Si > 0. 

Since Di is the support of Ui, the result follows. 



Ga2(x) 



□ 



Proof of Lemma 2.9 Let / be any continuous function which is superhar- 
monic on and > 7. By Lemma |2.2[ iii), s is harmonic on D, so f — s is 
superharmonic on D and attains its minimum in D on the boundary. Hence 
f — s attains its minimum in i7 at a point x where s{x) = 7(x). Since / > 7 
we conclude that f > s on 0, and hence everywhere. Thus s is at most the 
infimum in (22 ). Since the infimum in (22 ) is taken over a strictly larger set 
than that in (20), the reverse inequality is trivial. □ 
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7.2 Boundary Regularity 

For the proof below, we follow the sketch in Caffarelli [H Theorem 2]. 



Proof of Lemma 2.10. Fix xq G dD, and define 

L{x) = 7(xo) + (V7(a;o),x - xo). 



Since o" is C^, we have that 7 is by Lemma 2.3 ii). Let A be the maximum 
second partial of 7 in the ball B = i?(xo,4e). By the mean value theorem 
and Cauchy-Schwarz, for x €z B we have 

\L{x) -7(x)| = KV7(xo) - V7(x=^),x - xo)| 

< AVd\xo — — xol < ce^, (97) 

where c = 16 A\/^. Hence 

s{x) > 7(x) > L{x) — ce^, x G B. 

Thus the function 

w = s — L + ce^ 

is nonnegative and superharmonic in B. Write w = wq + wi, where wq is 
harmonic and equal to w on dB. Then since s{xo) = 7(xo), we have 

wo{xo) < w{xo) = s{xo) - L{xo) + ce^ = ce^. 

By the Harnack inequality, it follows that wq < c'e^ on the ball B' = 
B{xo,2e), for a suitable constant c'. 

Since wi is nonnegative and vanishes on dB, it attains its maximum 
in S at a point xi in the support of its Laplacian. Since Awi = As, by 
Lemma |2.2[ iii) we have s{xi) = j{xi), hence 

wi{xi) < w{xi) = s{xi) — L{xi) + ce^ < 2ce^, 



where in the last step we have used (97). We conclude that < w < 
(2c + c')e^ on B' and hence |s — L| < (c + c')e^ on B'. Thus on B' we have 

|u| = |s - 7I < |s - L| + I7 - L| < (2c + c')e^. 

In particular, u is differentiable at xq, and Vu(xo) = 0. Since s is harmonic 
in D and equal to 7 outside D, it follows that u is differentiable everywhere, 
and off dD. 
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Given y £ dD^^ let xq be the closest point in dD to y. Since -B(y, e) C D, 
by Lemma |2.2[ iii) the function w is harmonic on B{y, e), so by the Cauchy 
estimate [3,, Theorem 2.4] we have 

(J 

\Vs{y) - V7(xo)| = |V-u;(y)| < — sup w{z). 



Since B{y, e) C B' , the right side of (98) is at most C{2c + c')e, hence 

|Vn(y)| < |Vs(y) - V^{xq)\ + \V-i{y) - Vi{xo)\ < Coe 
where Co = C(c + c') + Thus u is on 9L> as well. □ 



We now turn to the proof of Proposition |2. 11 The first part of the proof 
follows Friedman ||10| Ch. 2, Theorem 3.5]. It uses the Lebesgue density 
theorem, as stated in the next lemma. 



Lemma 7.1. Let A C 



be a Lebesgue measurable set, and let 



A' = {xeA \ limmf e-'^C{B{x,e) n A^) > O}. 



Then C{A') = 0. 



To prove the first part of Proposition 2.11, given a boundary point x G 
dD, the idea is first to find a point y in the ball B{x, e) where u = s — j 
is relatively large, and then to argue using Lemma 2.10 that a ball B(y, ce) 
must be entirely contained in D. Taking A = dD in the Lebesgue density 
theorem, we obtain that x £ A'. 



The proof of the second part of Proposition 2.11 uses Green's theorem 
in the form 



D' 



{uAh — Auh)dx 



dD' 



dh du 
dn d 



n 



h dr. 



(99) 



Here D' is the union of boxes x^ that are contained in D, and dx is the 
volume measure in D' , while dr is the (d — 1) -dimensional surface measure 
on dD' . The partial derivatives on the right side of (99) are in the outward 
normal direction from dD' . 



Proof of Propostion 2.11. (i) Fix < A < 1. For x G d D wit h cr{x) < A, for 

-+2 

2 



small enough e we have a < on B{x, e). By Lemma 



2.5 ii) the function 



fiy) = i{y) + 



1 - A, 



\y\ 



64 



is superharmonic on B{x, e) n D, so by Lemma |2.2[iii) the function 



w{y) = u{y) 



1 - A, 



\x - y\ 



s{y)-f{y) + il-\){x,y) 



1 - A, 



is subharmonic on B{x, e)riD. Since w{x) = 0, its maximum is not attained 
on dD, so it must be attained on dB{x, e), so there is a point y with |x— y| = e 
and 



u{y) > 



1 - A 



By Lemma 2.10 we have |Vu| < (1 + c)Coe in the ball B{y,ce). Taking 
c = we obtain for z G B{y, ce) 

u{z) = u{y) + {Vu{y^),z - y) 

> ^^e' - c(l + c)Coe^ > 0. 

Thus for any x £ dD Ci {a < X} 

C{B{x, (1 + c)e) n {dDY) > ujaice^. 
By the Lebesgue density theorem it follows that 

C{dDn{a < X}) = 0. (100) 



By Lemma 2.7 we have o" < 1 on dD. Taking A | 1 in (100), we obtain 
£{dD) < £{dD n{a< 1}) + £(^-^(1)) = 0. 
(ii) Fix 5 > and let 

D'= U 



where x'-' = x + [— |, |]'^. By Lemmas 



2.2 iii) and 



2.3 



'ii), in D' we have 



Au = As - A7 = 2d{l - a) 



Now by Green's theorem (99 ), since u is nonnegative and h is superharmonic 

dh du 



1 — a)hdx = — 
o'' 2d 



Aun dx < — - 
~ 2d 
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where n denotes the unit outward normal vector to dD' . By Lemma 2.10 
the integral on the right side is bounded by 



dD' 



u 



dh 
dn 



+ 



du 
dn 



dr < CoVd{6''\\Vh\\^ + 6\\h\UCd-i{dD'), (102) 



where Cd-i denotes {d — l)-dimensional Lebesgue measure. Let 

Since C{dD) = 0, given e > we can choose 6 small enough so that C{S) < e. 
Since dD' C dS, we have 



2d£{S) 



Since D is open, C{D') | C{D) as (5 | 0. Taking 5 < e smaller if necessary 
so that C{D') > C{D) - e, we obtain from fm\ and ([Ml 

(1 - a)hdx < (M + l)e + Co\/d(||/i||oo + ||V/i||oo)e 



D 

where M is the maximum of |c7|. Since this holds for any e > 0, we conclude 
that h{x)dx < h{x)a{x)dx. □ 

Proof of Proposition \2.1S\ (i) Fix e > 0. Since a is continuous almost ev- 
erywhere, there exist functions (Tq < a < ai with /]gd(o'i — CTQjdx < e. 
Scaling by a factor of 1 + (5 for sufficiently small 6, we can ensure that 
C{a~^{l)) = 0, so that (Tq and ai satisfy the hypotheses of Proposition 2.11 
For i = 0, 1 let 

7j(x) = -|x|^ - Gai{x), 

and let Sj be the least superharmonic majorant of 7^. Choose a with A < 
a < 1 such that £{a^^{a)) = for i = 0, 1, and write 

Di = {si > 7i}; 

S^ = DiU{ai>a}. 



By (24) we have {uo > a} C {a > 1}°, hence by Lemma 2.8 
So = DoU {ao > a} C D U {a > 1}° = D 



(103) 



C D C DiU{a >1} = Si. (104) 
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For i = 0, 1 write 

ycFi[x), X f Di. 
By Proposition 2.11[ ii) with /i = 1, we have 



/ a° = £(A) + / a^= I 



For < ao < a, we have 

C{Si - So) < £{Si n {a°o < ao}) + C{{ao < < a}). 
Since uj' > a on Si, the first term is bounded by 

C{Sin{a°o < ao}) < 



(105) 



(106) 



Q'l ~ '^olli _ Iki ~ Q"o||i ^ e_ 



a — ao 



a — ao 



a — ao 



where the equaUty in the middle step follows from (105). 

By Proposition 2.11[ i) we have C{dDi) = 0, hence C{dSi) = 0. Taking 
ao = a — -y/e we obtain from (103), (104), and (106) 

C{db) <C(^i-So) =C{Si-So) < Ve + 'C({a-^/e < ao < a}). 

Since this holds for any e > 0, the result follows, 
(ii) Write 

^1, xeD 



a°{x) 



(t(x), X ^ D. 



From (105) we have 



/ ao= [ al< f a°< f = [ 

jRd J^d J^d J^d J^c 



CJl. 



The left and right side differ by at most e. Since J^a co ^ J^d o" < J^d ci, 
and e > is arbitrary, it follows that J^^ a = J^a o"° , hence 



[ a = I a° - ( a = C{D). 

Jd JRd Jdc 



(107) 



(iii) By Lemma 2.7, if (7{x) > 1 and x ^ D, then a is discont inuo us at x. 
Thus fj = 1 almost everywhere on D — D, and we obtain from (107) 



^a = C{D)+ a = C{D). □ 

D J D-D 
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7.3 Convergence of Obstacles, Majorants and Domains 

Proof of Lemma 2.16\ (i) Since a, a-n are supported on we have for x G 



Gan{x) - G(j{x)\ < I \g{x,y)\\an{y) - cF{y)\dy. 



(108) 



Fix < e < 1. Since < a, (t„ < M, we have by (|T7|) and (|T8 
/ \9{x,y)\\an[y) - (T[y)\dy < I 

JB{x,e) I 



j2Me2 (21og^ + l) , d = 2 



2d 
[. d-2 



d>3. 



(109) 



Now let Bq = B — B{x, e). Since (T„ — > a in L^, we can take n large enough 
so that |o"n — 0"! < e"'"^, obtaining 



|5(x,y)||cj„(y) -(7(y)| < e'^ ^uj^R'^ sup 15(2;, y)| 
Bo yeBo 



< 



'2i?2emax(logi,log(i?+ |x|)), d = 2; 



d=2R ^' 



d > 3. 



Thus the right side of (108) can be made arbitrarily small for sufficiently 
large n. 

(ii) By Lemmas 



2.8 



and 



2.14 



the non coincidence sets D, are con- 



tained in the ball Bi = M^^'^B. Given e > and a compact set K containing 
Bi, choose A'" large enough so that |7„ — ^\ < e on K for all n> N . Then 

s + e>7 + e>7rt 

on K, so the function fn = max(s + e, 7^) is superharmonic on K. By 
Lemma 2.9 we have fn > s„, and hence s + e > Sn on K . Likewise 

Sn + e>7n + e>7 

on s o th e function /„ = max(s„ + e, 7) is superharmonic on K. By 
Lemma 



2.9 



we have fn^s and hence + e ^ s on ii'. Thus s„ — > s 
uniformly on K. 

(iii) Let /? > be the minimum value of s — 7 on D^, and choose large 
enough so that I7 — 7„|, \s — Sn\ < /3/2 on D for all n > N. Then 

Sn-7n>s-7-/?>0 



on De, hence C ^^(n) for all n > 



□ 



68 



Proof of Lemma 2.17, Consider the functions 

/±(2/)=/(y)±A|y|2. 



(110) 



By hypothesis, /+ is subharmonic and /_ is superharmonic on B{o,R). 
Given x £ B(o,f3R), let r = min(i?, 2|x|), and let h± be the harmonic 
function on the ball B = B(o,r — 1) which agrees with f± on dB. Then 



f+<h+<h^ + 2Ar2 < /- + 2Ar^ 



(111) 



By the Harnack inequality |19| Theorem 1.7.2] there is a constant such 
that 

h+{x) < Cf3h+{o). 



Since f~{o) = f{o) = we have from (111 ) 

h+{o) < f-{o) + 2Ar2 < 8A|xp. 

hence 



f{x) < /i+(x) — A|a;| < c^A|x| 



with c/3 



1. 



□ 



Lemma |2.18 is proved in the same way as Lemma 2.17 using the con- 
tinuous Harnack inequality in place of the discrete one, and replacing |yp 
by \y\'^/2d in (fTTol). 



Proof of Lemma 2.19. Let K be the closure of (^Ui3)e, and let Yn = >1]^/„U 
Bi/n- Since K is contained in |J 1^ = yl U -B, the sets K CiYn form an open 
cover of K, which has a finite subcover, i.e. K C Yn for some n. □ 



7.4 Discrete Potential Theory 



Proof of Lemma 2.20. By Taylor's theorem with remainder 



d"^ f I d^f 

f{x+6nei)-2f{x)+f{x-6nei) = ^"^n + g^' 



6 dx^ 



{x-uei)6f, 



for some < t, u < 1. Summing over i = 1, . . . , d and dividing by 6^ gives 
the result. □ 
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Proof of Lemma 2.21. In dimensions d > 3 we have from (25) and the stan- 



dard estimate for the discrete Green's function [30] (see also |19| Theorem 
1.5.4]) 



gn{x,y) =Sl lad 



\x - y\ 



2-c 



+ 



\x - y\ 



= ad\x-y\^-'' + 0{6l\x-y\-''). 
Likewise, in dimension two, using the standard estimate for the potential 



kernel ^T] (see also |19| Theorem 1.6.2]) in (26) gives 



, > 2 ^ \x — y\ 2 , . ^ f\x — y\ 
9n{x, y) = — log — h - log (3„ + O — - — 

TT On Vr \ On 



□ 



Proof of Lemma \2.2^ Let C M*^ be compact. By the triangle inequality 
|(G„cT„)° - Ga\ < \{Gnan)° - Ga°\ + \Ga° - Ga\. 



By Lemma |2.16[ i) the second term on the right side is < e/2 on X for 
sufficiently large n. The first term is at most 



\{Gnan)°{x)-Ga°{x)\< ^ cT„(y) 



^n9n{x-',y) - / g{x,z)dz 

a 



^ i \9ni.x-",y) - g{x, 



z)\dz 



M I \gn{x-'- , z'-'-) - g{x,z)\dz. 



(112) 



By Lemma |2.21[ we have 



\gnix'-'- , z'-'-) - g{x, z)\ < \gn{x-'- , z'-'-) - g{x-'- , z-'-)\ + \g{x'-'- , z'-') - g{x, z)\ 
< C6Jx-z\'^-'^ 



for a constant C depending only on d. Integrating (112) in spherical shells 
about X, we obtain 



|(G„a„)° - Ga°\ < Cduj^MRSn, 

where R is the radius of B. Taking n large enough so that the right side is 
< e/2, the proof is complete. □ 
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Figure 4: Modified rotor-router aggregation started from a point source at 
the origin in Z^. In the left figure, particles are killed on hitting the positive 
X-axis. In the right figure, sites on the positive x-axis send all particles in 
the downward direction. 



Proof of Lemma 2.23. Let / be any function which is superharmonic on i7 
and > 7„. Since Sn is harmonic on D„, the difference f — Sn '^s superharmonic 
on Dn and attains its minimum in Dn U dDn on the boundary. Hence f — Sn 
attains its minimum in 17 U dfl at a point x where s„(x) = 7n(x). Since 
f ^In'^^ conclude that / > on Q and hence everywhere. Thus s„ is at 
most the infimum in ( 38 ) . Since the infimum in ( 38 ) is taken over a strictly 
larger set than that in (36), the reverse inequality is trivial. □ 



8 Open Problems 

We conclude by mentioning two simple variants of aggregation models which 
appear to have interesting limit shapes, but about which we cannot prove 
anything at present. In the first variant, particles start at the origin and 
perform rotor-router walks in J? until either reaching an unoccupied site or 
hitting the positive x-axis. The resulting set of occupied sites is shown on 
the left in Figure |4j In the second variant, shown on the right, particles 
again start at the origin and perform rotor-router walks until reaching an 
unoccupied site, except that the rotors on the positive x-axis do not rotate 
but instead always point downward. 

Our simulations indicate that limiting shapes exist for both of these 
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variants. Pavel Etingof (personal communication) has found exact solutions 
to the analogous Hele-Shaw flow problems in the continuum. The challenge 
is to extend our methods to prove that the lattice models produce shapes 
converging to these Hele-Shaw solutions. 
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